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Abstract 

We consider the effective theory of large D stationary black holes. By solving the Einstein 
equations with a cosmological constant using the 1/D expansion in near zone of the black 
hole we obtain the effective equation for the stationary black hole. The effective equation 
describes the Myers-Perry black hole, bumpy black holes and, possibly, the black ring solu¬ 
tion as its solutions. In this effective theory the black hole is represented as an embedded 
membrane in the background, e.g., Minkowski or Anti-de Sitter spacetime and its mean 
curvature is given by the surface gravity redshifted by the background gravitational field 
and the local Lorentz boost. The local Lorentz boost property of the effective equation 
is observed also in the metric itself. In fact we show that the leading order metric of 
the Einstein equation in the 1/D expansion is generically regarded as a Lorentz boosted 
Schwarzschild black hole. We apply this Lorentz boost property of the stationary black 
hole solution to solve perturbation equations. As a result we obtain an analytic formula 
for quasinormal modes of the singly rotating Myers-Perry black hole in the 1/D expansion. 
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1 Introduction 


The infinite dimensional limit of General Relativity gives not only various interesting 
pictures of black hole physics, but also new useful analytic method to solve gravitational 
problem [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. The perturbation problem, such as quasinormal 
modes (QNMs), has been solved for general static black holes [7, 8] and the Myers-Perry 
black hole with equal spin [6]. As for non-linear problem the effective theory for the large 
D black hole has been considered recently in [9, 10] and some non-trivial static solutions 
were constructed. In this paper we study the effective theory of the large D stationary 
black hole in asymptotically flat or Anti-de Sitter (AdS) background. 

The reason for the simplification of the gravity at large D is simple [2]. The radial 
gradient of the gravitational field around black holes becomes very large with 0{D/tq) at 
large D. ro is the black hole size. So the tangential dynamics along the horizon becomes 
sub-dominant, and the system reduces to the ordinarily differential equation system with 
respect to the radial direction. In addition the appearance of the hierarchical sale, vq/D 
and ro, gives two separated excitations around the black hole, that is, the non-decoupled 
mode with the frequency euro = 0(D) and decoupled mode cor q = O(D 0 ). While the non- 
decoupled mode is universal for black holes [5], the decoupled mode has a particular and 
interesting feature for each black hole. For instance the instability of the ultraspinning 
back holes and black branes are belonging to this decoupled sector [2, 6]. These two facts, 
the dominance of the radial dynamics and appearance of two hierarchical scales, make 
General Relativity simple and analytically solvable, but still non-trivial at large D. 

Achievements in QNMs by large D method motivate us to construct the effective theory 
for general black hole solutions in the 1/D expansion as a next step. The gravitational 
field of the black hole is localized in the near region of the horizon due to its large gradient 
at D —>• oo. Then the black hole can be regarded as a membrane in the background 
spacetime like the membrane paradigm [11, 12, 13] at this limit. The physical feature 
of is determined by solving Einstein equations in near zone of the black hole. If we 
consider the decoupled mode dynamics of the black hole, the Einstein equations can be 
solved consistently only in near zone. The obtained solution gives the boundary condition 
for Tib as physical properties. Then we can construct the effective theory of the large 
D black hole as the membrane E^ in the background spacetime. This construction of 
the effective theory has been considered in [9] for static black holes and for more general 
solutions with time dependences in [10]. 

Here we study the effective theory for general large D stationary black holes in the 
Minkowski or AdS background such as the (AdS)Myers-Perry black hole, bumpy black 
holes and black ring. These solutions have been known as exact solutions [14, 15, 16] or 
numerical solutions [17, 18, 19]. We obtain the general analytic metric for the stationary 
black holes in the 1/D expansion, which describes these solutions in the unified manner. 
The solution gives the equation for the embedding of E^ in the background. This embed- 
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ding determines the horizon topology and its geometric shape. Then the effective equation 
for the embedding of 'Eb is given by the equation for the mean curvature K of as 


7 


-l 


K 


Ss 


= 2k, 


( 1 . 1 ) 


where 7 is a redshift factor containing both the gravitational redshift on the membrane 
from the background geometry and the Lorentz redshift from the local rotation along the 
membrane, i.e., 


7- 1 = (1.2) 

k and Qh are the surface gravity and horizon angular velocity of the black hole respectively, 
is the rotation radius of Eg defined by 


T^<j) = 



(1.3) 


A solution of this effective equation gives the black hole metric which solves the leading 
order Einstein equation in the 1/D expansion. This effective equation is one of main 
results of this paper. In the static case this effective equation reduces to one of static 
black holes [9]. 

The effective equation implies that the mean curvature of E# is regarded as the surface 
gravity redshifted by the local Lorentz boost. This local Lorentz boost property of the 
effective equation turns out to be also the property of the leading order metric itself. 
Actually we find that the leading order metric in near zone is obtained as a Lorentz 
boosted Schwarzschild black hole. This property of the leading order solution would not 
hold if we go to the higher order structure in the 1/D expansion. However it can still 
play a crucial role in solving the Einstein equations. As one interesting example, we 
show that perturbation equations on the obtained leading order general stationary black 
hole can be explicitly solved by using the solution of one of the Schwarzschild black hole. 
This perturbation solution, of course, contains the solution of perturbation equations of 
the Myers-Perry black hole. Going to the higher order of perturbation equations in the 
1/D expansion we obtain an analytic formula for QNMs frequency of the singly rotating 
Myers-Perry black hole both for the axisymmetric and non-axisymmetric perturbations. 
As a result the threshold angular momentum of the ultraspinning instability of the singly 
rotating Myers-Perry black hole is found to be 

a 2 

- = 2k s -l, k s = 1,2,3,... (1.4) 

r o 

where ro is the horizon radius and ks is a positive integer. This result can be obtained 
in a much easier and simpler manner by examining the effective equation ( 1 . 1 ) for the 
embedding of the bumpy black hole. Note that there is no dynamical instability at the 
threshold angular momentum of the case ks = 1. The threshold for ks = 1 does not imply 
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any existence of a bifurcation zero-mode. There is only a trivial perturbation along the 
Myers-Perry solution branch. 

The remaining of this paper is organized as follows. In section 2 we consider the 
effective theory for stationary black holes at large D. We give the leading order solution 
and its boost representation. The effective equation for the large D stationary black holes 
is derived from this leading order metric. We study a simple solution and non-trivial 
solutions of the effective equation there. In section 3 we apply the boost property of the 
leading order solution to solve perturbation equations on the general stationary solutions. 
The section 4 pursues the consideration for the perturbation in the previous section in more 
detail. We concentrate on the perturbation of the Myers-Perry black hole and obtain the 
QNMs frequency analytically. We close this paper by discussion and outlook in section 5. 
The appendices contain some technical details of calculations for the main part. 


2 Large D stationary black holes 


We study the D dimensional stationary black hole solution and its large D effective theory. 
The following analysis is performed in the similar way with static case [9] 1 . We use the 
small expansion parameter 1/n where 


instead of 1/D. 


n = D — 3, 


( 2 . 1 ) 


2.1 Set up 


The large D black hole has a very large gradient along the radial direction, p, with O(n). 
Then our metric ansatz for the D dimensional stationary black hole is given by 

ds 2 = dp 2 _|_ g ab dx a dx b . (2.2) 

n z 

z is one of D — 1 dimensional coordinate x a . In this paper we consider the metric which 
has only one inhomogeneous coordinate, z. other than p for simplicity. The generalization 
to several inhomogeneous coordinate case would be straightforward as done in [9]. Then 
the Einstein equation with a cosmological constant 


(D-l)(D-2) 

L 2 


is decomposed on p = constant surface to 


—R + K 2 — K a h K b a 
V a K a b - V b K = 0, 


(D-l)(D-2) 

L 2 


= 0 , 


71 JD — 1 

-d p K% = -KK\ + R\ + -p-n - - V a V b N, 


(2.3) 


(2.4) 

(2.5) 

( 2 . 6 ) 


But the index notation is changed from [9] 


i 
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and 


K °- = (2 - 7) 

R a b and V a is the Ricci curvature and covariant derivative of D — 1 dimensional metric 
g a b respectively. K a b is the extrinsic curvature of D — 1 dimensional p = constant surface. 

We consider stationary solutions with a rotation. Then the D — 1 dimensional metric 
ansatz for g a b is 

= g a bdx a dx b 

= —A(p, z) 2 dt 2 + F(p, z) 2 (dcj) — W(p, z)dt) 2 

+ G(p, z) 2 dz 2 + H(p, z) 2 qABdx A dx B , (2.8) 

where x A and qab are a coordinate and the standard metric on S D ~ A respectively. We 
solve the Einstein equations for K a b and metric functions A(p, z ), F(p , z), W ( p , z), G(p, z) 
and F[(p, z) by using the 1/D expansion. In the following we specify large D behaviors of 
these functions and boundary conditions as a set up. 

Large D behavior At first each metric functions should be O(n 0 ). In particular 

g zz = O(n 0 ) (2.9) 

implies that the derivative with respect to z is O(n 0 ). The hierarchical scaling between p— 
and £— dependence makes the leading order equation at D —>• oo be ordinary differential 
equations with respect to p. Next we assume that the leading order of G(p, z) and H(p, z) 
has only ^-dependence and no p-dependence. This assumption is motivated by the exact 
solution of the Myers-Perry black hole (see Appendix A) and the fact that the large 
dimensionality of the sphere gives another D factor in p-dependent terms through the 
trace on S D ~ A . Thus we assume 

K% K a b = 0(n°). (2.10) 

while the leading order of A(p,z), F(p,z) and W(p,z ) has p-dependences as 

K t u K t ^K^K = 0(n). (2.11) 

From these assumptions we can derive the scaling law of various geometric quantities at 
large D. For example, the Ricci scalar R of g a b is 0(n 2 ). 

Boundary condition The boundary conditions are imposed on the horizon and at the 
asymptotic infinity of near zone defined by n > p > 1 At the horizon the extrin¬ 
sic curvature has a singular behavior essentially only in one component. This condition 

2 The effective theory of the large D black holes is described as the membrane physics seen from the 
far zone [9]. In the overlap zone n > p > 1 we can obtain physical properties of the membrane by the 
near zone solution as the boundary condition. 
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guarantees the regularity of following quantities 

K-K\-K\, K* z , K a b , (2.12) 

on the horizon. For the static case K''l should be regular on the horizon. In this paper 
we consider the asymptotically flat or AdS black hole. Then the boundary condition at 
asymptotic infinity of near zone is 

A(p,z) = \/ 0 (z) + 0(l/n,e~ 2p ), W{p, z) = 0(e~ 2p ). (2.13) 

Note that the Newton potential of the black hole is 0(e~ 2p ) in our gauge. So this asymp¬ 
totic boundary condition means the leading order metric becomes the flat or AdS metric 
other than the potential term at p> 1. While g m should be damping exponentially in p 
at all order of 1 /D, gtt has some terms which does not decay at the asymptotic infinity 
in the higher order solution in the 1/D expansion when we consider the asymptotically 
AdS case (See Appendix A). Since we solve the leading order equations in this section, we 
do not explore this structure in detail. The boundary condition (2.13) is enough for the 
leading order solutions. 


2.2 Solving the leading order equation 

We solve the leading order equation of the Einstein equation in the 1 /D expansion by per¬ 
forming p integrations. This integration gives free functions of 2 as integration functions. 
They consist of the effective theory of large D stationary black holes. 

Our large D scaling assumptions give following large D behaviors 


N(p, z) = Nq(z) + 0(l/n), (2.14) 

G(p,z) = l + 0(l/n), (2.15) 

H(p,z) = TZu(z) + 0(l/n), (2.16) 


where p-dependence of the leading order of N(p, z) and z-dependence of the leading order 
of G(p, z ) have been absorbed into the definition of p and z coordinate respectively. The 
leading order of the Ricci scalar of g a b can be calculated from this ansatz as 


„ z . 1 i n'izf 
~ n 1 1? + ~ 


(2.17) 


where here and in the following we omit 0(l/n) symbols showing the existence of sub¬ 
leading corrections in 1 /D for the simplicity in the representation. Then the trace part of 
eq. (2.6) is integrated to 


K = T"T c °th (^Tjip - Po{z))\ 
ro(z) V r o(z) J 

where we defined tq{z) by 

1 1 . 1 no*) 2 

r 0 {zY Li + K n (z ) 2 K k {z) 2 ' 


(2.18) 


(2.19) 
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Po(z) is the horizon position p = po(z), so we can set to po(z) = 0 by a gauge choice. The 
residual gauge in p coordinate can be used to have 


N 0 (z) = 2 r 0 (z). 


( 2 . 20 ) 


Note that this choice for Nq(z) is different from [9] by factor two. This is just a gauge 
choice and it does not affect physical properties. The equations for K t t , and K° 0 in 
eq. (2.6) are solved by 


K t = Ctt(z) t = Ct<j>(z) $ = C<t><p(z) 

* ro(z) sinh 2p’ ^ r 0 (z) sinh2p’ ^ ro(z) sinh2p’ 


( 2 . 21 ) 


Ctt(z), Ct(j,(z ) and C^^z) are integration functions with respect to the p integration. The 
boundary condition on the horizon (2.12) gives 


1 - Cu(z) - C^j,{z) = 0. 


( 2 . 22 ) 


In the following we eliminate C^Az) by this condition. Then we can integrate other 
components of eqs. (2.6) and (2.7), and obtain the following leading order solutions 


and 


A(p, z) 2 
F(p,z) 2 
W(p,z) 


Aq(z) 2 Fq(z) 2 tanh 2 p 
Fq(z) 2 - Ao(z) 2 Ctrf,(z) 2 tanh 2 p’ 

F 0 (z) 2 - A 0 (z)C t Az) 2 tanh 2 p, 

Tb(;Q 2 (l - C u (z)) - A 0 (z) 2 C t Az) 2 C a (z) tanh 2 p 
C t <f,(z)F 0 (z) 2 - A 0 (z) 2 C t Az) 3 taiAi 2 p 


(2.23) 

(2.24) 

(2.25) 


G(p,z) = 1 - 


2ro(z) 2 (IZ'Az) 1 \ 

n S( P) ’ 


H(p,z) = 7l H (z) l 1 + - log (cosh p) 

n 


(2.26) 

(2.27) 


Aq(z) and Fq(z) are integration functions. The integration functions of G(p, z ) and F[(p, z) 
are absorbed into the 0{\/n) redefinition of z coordinate and 77 h(z)- The scalar constraint 
equation (2.4) could be satisfied by choosing the integration function in W(p, z) properly. 
We can see that these solutions reduce to the static solution [9] if we set to 


Ctt(z) — 1, Cty — 0. 


(2.28) 


We could solve the scalar constraint equation (2.4) and all evolution equations of eqs. 
(2.6) and (2.7). The remaining equation is the vector constraint equation (2.5). The 
vector constraint equation gives an additional condition between integration functions as 



( A 0 (z)C tt (z) \ 
V r 0 (z) ) 


(1 - C tt (z))— log 


( Ctt(z) 

\Ct4>(z) 


= 0 . 


(2.29) 


This equation, at first glance, seems to be a non-trivial condition. However we can make 
this equation trivial in the following sense. To realize this we observe the surface gravity 
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k and angular velocity fl// of the horizon of the leading order solution. These quantities 
are read as 


_ nAo(z) 
K> — ~ ■ 


Qh = 


1 - c tt (z) 


(2.30) 


2 r 0 (z) ’ C t(j> {z) 

These quantities should be constant on the horizon [20]. Then, using the condition that k 
and Qfj are constant, we can see that the constraint equation (2.29) is satisfied automati¬ 
cally. Thus the vector constraint equation becomes trivial under the constancy condition 
of the surface gravity and angular velocity of the horizon. This is the same situation 
with the static case [9]. In the static case the vector constraint is equivalent to the con¬ 
stancy condition of the surface gravity. In stationary case the vector constraint can be 
satisfied if we use the constancy condition of not only the surface gravity but also the 
horizon angular velocity. One may feel that this statement is strange since originally the 
constancy of the surface gravity and angular velocity of the horizon was shown by using 
the Einstein equations [20]. It means that the constancy condition of the surface gravity 
and horizon angular velocity should be derived conditions, not additionally imposed one 
to satisfy the Einstein equation. In fact the constancy conditions (2.30) can be derived 
if we use eq. (2.29) and another equation which is obtained at the next-to-leading order 
of 1 /D expansions. The regularity condition of the next-to-leading order correction of 
K z z on the horizon requires one additional condition between Aq(z ) and Fq(z) (see eq. 
(B.17)). Then these two equations give the constancy condition of the surface gravity 
and horizon angular velocity as eq. (2.30). The detail of this argument will be given in 
the appendix B. So, here, we assume the constancy conditions (2.30) on the leading order 
solution in advance. If constancy conditions (2.30) are not satisfied, the next-to-leading 
order solutions becomes singular at the horizon although the leading order solutions still 
regular. So we can replace Aq(z) and Ct<f>(z) by Ctt{z), k and using eq. (2.30). 

Finally we impose the boundary condition at the asymptotic infinity (2.13) on the 
leading order solutions. Then we find 


C tt (z) = 


V 0 (z) 2 „ , ,2 V 0 (z) 2 — 4/c 2 r 0 (z) 2 


r, F 0 (zY = 


4« 2 ro(z) 2 ’ uv ~ y Ak 2 Q 2 H ro(z) 2 

where we defined the reduced surface gravity k by 


k = nn. 


As a result we obtain regular leading order solutions as 

4k 2 ro(z) 2 \/o(z) 2 sinh 2 p 


A(p,z) = 


Vq (z) 2 + 4k 2 ro(z) 2 sinh 2 p' 


F(p,z) 2 = (y 0 (z) 2 -Ak 2 r 0 (z) 2 ) 


2 . Vq(z) 2 + 4:k 2 ro(z) 2 sinh 2 p 


4Fl? H k 2 ro{z) 2 cosh 2 p 


W(p,z) = 




Vq (z) 2 + 4k 2 ro(z) 2 sinh 2 p' 


(2.31) 

(2.32) 

(2.33) 

(2.34) 

(2.35) 
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and eqs. (2.26) and (2.27) for G(p,z ) and H(p,z). One can check that the (AdS)Myers- 
Perry black hole is described by this leading order solution as seen in Appendix A. 


2.3 Effective equation 


In the same spirit with [9] we can obtain the effective equation for the stationary black 
hole from the leading order solution of near zone obtained above. The near zone metric 
is matched with the far zone metric at overlap zone n 3> p 1 under the boundary 
condition (2.13). The near zone metric on p = constant surface T,b at overlap zone is 


j _2 i \/ / _.\2 , Vo(z) 2 - 4:k 2 r 0 (z) 2 

ds Is# vo(^) dt + 2 dcj> 

\i H 


+ dz 2 + TZb{zY qABdx A dx B , 


(2.36) 


where we neglect the potential terms 0(e 2p ). Notice that the trace of the extrinsic 
curvature, K , is read as 


K 


Ss 


n 


(2.37) 


where ro(z) is given in eq. (2.19). The gravitational redshift factor of by the back¬ 
ground spacetime, gtt , is 


\f—Qtt |s fl = V 0 (z). 


(2.38) 


Using eqs. (2.37) and (2.38), we can see that the leading order solution in near zone gives 
the following embedding equation for E# at overlap zone as 




= 2k, 


where we defined the rotational radius by 


(2.39) 


TZ(j) = 



(2.40) 


The effective equation (2.39) gives the mean curvature of as the surface gravity red- 
shifted by the gravitational effect and the local Lorentz boost effect with the boost velocity 
defined by 


Vcj) = (2-41) 

The regular leading order solution contains arbitrary functions, Vo (z) and IZn(z). These 
functions can be determined by specifying the embedding into the background. If one 
consider the embedding into the Minkowski background, we have Vo(z) = 1 and there is 
only one remaining function IZu(z) which should satisfy eq. (2.39). On the other hand, if 
one considers the embedding into AdS background, the embedding would give a non-trivial 










redshift factor Vq (z) of £#. Thus the embedding into the AdS background should satisfy 
eqs. (2.38) and (2.39) for Vo (z) and 1Zy{(z). 

The effective equation (2.39) for the stationary black hole can be still interpreted by 
a soap-film equation as one for the static case [9] 3 . The membrane (soap-film) between 
two fluids satisfies the following equation for the mean curvature K of the membrane (see 
eq. (2.15) in [21]) 


K = 


yTAs — A p 
a 


(2.42) 


where a is the surface tension, T is the temperature, As and A p are the differences between 
the entropy density and the energy density of the fluids at both sides of the membrane. 
7 is the same redshift factor as in eq. (1.2): in eq. (2.14) in [21], one just has to extract 
£ 2 = —gtt- If we now consider a case in which A p is negligible, then the equation (2.42) 
becomes 


7 ~ X K =-, (2.43) 

<7 

so if we identify 2k = T As/a, then this is the same as eq. (2.39) (of course in this iden¬ 
tification T, As, and a are not independently determined, instead only their combination 
into k is fixed). 

There are two remarks on the metric (2.36) at £#. One is that the metric (2.36) does 
not need to have a component of gt The metric (2.36) describes the background metric 
where the solution is embedded, and it does not need to be rotating itself although the 
solution is rotating. Thus the metric (2.36) does not need to have gt<p component. Another 
remark is about the static limit. The static limit corresponds to the zero horizon angular 
velocity Qh = 0. Then the metric (2.36) requires 


V 0 (^) 2 = 4 k 2 r 0 (z) 2 


(2.44) 


at the static limit. One may think that the appearance of this additional condition is 
curious. However one can immediately find that the condition (2.44) is nothing but the 
effective equation of general static solutions [9] 


sTmK 


Sfl 


2k, 


(2.45) 


by using eqs. (2.37) and (2.38). This embedding equation can be obtained also by the 
static limit of eq. (2.39). So there is no appearance of new additional condition at the 
static limit. 

The embedding in the Minkowski background has a constant redshift factor as 


Vo(z) = 1. 


(2.46) 


3 This interpretation of the effective equation (2.39) is due to Roberto Emparan. We thank him for 
sharing it with us. 
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For the static embedding in this Minkowski background we have only one unique embed¬ 
ding by a round sphere corresponding to the Schwarzschild black hole [9, 22 ], This fact 
is consistent with the uniqueness theorem of the static solution [23]. On the other hand, 
however, the embedding of the stationary solution allows various embeddings as seen be¬ 
low. This is due to the appearance of new degree of freedom in the horizon deformation 
described by fin 77-0 in the effective equation. 

2.3.1 Ellipsoidal embedding 

As one application of the effective equation (2.39) we consider an ellipsoidal embedding 
in the Minkowski background. We embed the leading order solution by r = r{9) in the 
Minkowski background using the ellipsoidal coordinate. The flat metric in the ellipsoidal 
coordinate is 

u 2 j,2 r 2 + a? COS 2 9 2 , / 2 , 2 2 a \j a 2 

as = —at H-- 77 - dr + (r + a cos 9)d0 

r A a i 

+ (r 2 + a 2 ) sin 2 9d<p 2 + r 2 cos 2 9dVL 2 D _ 4 . (2.47) 

a is a parameter describing the oblateness of the ellipsoidal. We set 

« = (2-48) 

by using the normalization of near zone f-coordinate for simplicity. The relation z and 8 
can be seen from eqs. (2.36) and (2.47) with r = r(9 ) as 

dz I (1Zn(9) 2 + a 2 cos 4 9)(lZn(8) 2 + (1Z' H (9) + tan 9TZu(9)) 2 + a 2 cos 2 9) 
d9 \j cos 2 9(TZu(9) 2 + a 2 cos 2 9) 

where we used the relation 

r{8) cos 9 = K n {8) (2.50) 

in the leading order solution. Then the effective equation (2.39) becomes 

Q 2 h tan 2 6(a 2 cos 4 9 + TZu(9) 2 + cos 9 sin 9TZ}i(9)7l' }:{ (9)) 2 (a 2 cos 2 8 + 77h(0) 2 ) 

+ (77 h ( 0) 2 - l)(a 2 cos 4 9 + 1Zh( 9) 2 ) 2 
+ 2cos9sineiZn{8)n' n {9)(JZn{9) 2 - l)(a 2 cos 4 8 + 77 h (#) 2 ) 

+ cos 2 8n n {.8) 2 n' n {8) 2 {a 2 cos 4 8 + 77 h (#) 2 - sin 2 9) = 0 . (2.51) 

We solve this equation for TZn(9) in some cases below. 

Myers-Perry solution The Myers-Perry solution is described by the embedding 

7^h( 0) = cos 8 , Q H = ——- (2.52) 

1 + a z 

in the ellipsoidal embedding for arbitrary a. 
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Bumpy black holes The Myers-Perry black hole has marginally stable quasinormal 
modes u = 0 in the ultraspinning region [24, 25, 26]. The marginally stable modes suggest 
the existences of new deformed solution branching off the Myers-Perry solution branch. 
These solutions, so called bumpy black holes, have been constructed numerically in [17, 18]. 
Here we consider this solution by using the effective equation (2.51). To study this we 
perform the perturbative analysis of the effective equation (2.51) around the Myers-Perry 
black hole. We expand the embedding function 77 h( 0) around the Myers-Perry black hole 
as 


where 


77 h ( 0) = cos e\l + h{6)e + O (e 2 ) 


(2.53) 


e = tin 


a 

1 + a 2 


(2.54) 


Then, perturbing eq. (2.51) with respect to e, we find the perturbative solution 


n{e) 


a(l + a 2 ) 2 sin 2 6 sin 1+ “ 2 6 

(1 — a 2 )(l + a 2 cos 2 6) l + a 2 cos 2 0’ 


(2.55) 


where A is an integration constant describing the perturbation amplitude. The first part 
describes just the trivial perturbation adding the angular momentum along the Myers- 
Perry black hole solution branch. The second part is the non-trivial deformation of the 
solution into the bumpy black hole solution branch. For general a and 4 / 0, we have 
a non-analytic behavior in the solution at 9 = 0. The regularity of the solution at 6 = 0 
requires 4 


a 2 = 2/c — 1, As = 1,2,3,... 


(2.56) 


For k = 1, i.e., a c = 1, the solution is same with the first part of eq. (2.55), so it 
corresponds to the Myers-Perry black hole with a different angular momentum. The value 
a c for k = 1 is known analytically in all D by the thermodynamic argument [24] as 


D -3 


D -3 (D -3 


2(D — 4) \D — 5 


(D- 5)/2 


(2.57) 


ro is the horizon radius of the Myers-Perry black hole (see Appendix A). At the large D 
limit the formula (2.57) gives same value with eq. (2.56) for k = 1. For k > 1 we have non¬ 
trivial bumpy black hole solutions by eq. (2.55). Eq. (2.56) gives the threshold angular 
momentum of the ultraspinning instability of the Myers-Perry black hole. In Table 1 we 

4 Here the regularity means that any derivatives of TZ(6) should be finite at 9 = 0 and 7 Z(9) be an even 
function around 9 = 0. To satisfy this requirement a should be discretized by k as given in eq. (2.56). As 
we can see in appendix C this regularity condition at 8 = 0 is equivalent to the condition for the spheroidal 
harmonics at large D. Furthermore this discretization parameter k corresponds to the angular momentum 
number l. 
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Q j c\k=l 

&c | k=2 

CO 

II 

O 

<3 

&c | k =4 

eq. (2.56) 

l 

V3 = 1.732 

y/5 = 2.236 

V7 = 2.646 

numerical (D = 6) 

1.097 

1.572 

1.849 

2.036 

numerical (D = 7) 

1.075 

1.714 

2.141 

2.487 

numerical (D = 8) 

1.061 

1.770 

2.275 

2.725 

numerical (D = 9) 

1.051 

1.792 

2.337 

2.807 

numerical (D = 10) 

1.042 

1.795 

2.361 

2.855 

numerical (D = 11) 

1.035 

1.798 

2.373 

2.879 


Table 1: Comparison of our analytic formula (2.56) for the threshold of the instability 
of the Myers-Perry black hole with numerical results in [26]. The analytic formula shows 
quite good agreements with numerical results within the expected error 0(1/ D). The 
value of a c for k = 1 by [26] has perfect agreements with eq. (2.57). 


show the comparison of eq. (2.56) with numerical results in [26]. Actually we have good 
agreements within the expected error 0(1/0). This threshold angular momentum will be 
reproduced also by observing the quasinormal mode frequency directly in section 4. 

In Appendix A we generalize the above analysis to the AdS Myers-Perry black hole. 
Then we obtain the threshold angular momentum for the ultraspinning instability of the 
AdS Myers-Perry black hole in the 1/0 expansion. 


2.4 Boost representation 


One important property of the leading order solution in near zone is its local Lorentz boost 
representation. The effective equation (2.39) has the interpretation that the embedded 
surface Sb is the locally Lorentz boosted membrane. This local Lorentz boost property 
holds also for the metric itself in near zone. To see this we rewrite the (t, <f>) part of the 
obtained leading order metric in near zone as 


ds {t,ct>) = -Vo(^) dt + 


+ 


2 r .2 , V o(^) 2 -4K 2 r 0 (2) 2 2 


O 2 


4k 2 ro(z) 2 cosh 2 p 


-d,r 


Mo{z) 2 dt — 


Vo (z) 2 - 4/t 2 r l 

Q h 


d<f> 


(2.58) 


By introducing the boost parameter cr(z) defined by 


and new local frame dx p by 


cosh a(z) 


Vq(z) 

2kr 0 (z) 


dx p (z ) 


(\/ 0 {z)dt, 


y/V 0 (z) 2 - 4 k 2 r 0 (z) 2 



(2.59) 


(2.60) 
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the leading order metric is written in a very simple form as 

ds 2 = dp 2 + (r) pq + UpUq ) dx p dx q + dz 2 + TZu(z) 2 dQ 2 D _ 4 , (2.61) 

n z \ ™ cosh p) 

where r] pq is the flat metric on two dimensional spacetime. The fluid velocity u p is 

u p dx p = Vo (z) cosh a(z)dt — ^ ° i_J_ sinh a(z)d(/>. (2.62) 

Note that the leading order metric of the Schwarzschild black hole is (see Appendix A) 

ds 2 = — dp 2 — tanh 2 p dt 2 + dz 2 + sin 2 9d(p 2 + cos 2 9d£i 2 D _ 4 . (2.63) 

Thus the stationary solution can be represented as the locally boosted Schwarzschild black 
hole under the following boost transformation 

dt -> V 0 (z)dt cosh a(z) - V v °(^) 2 ~ 4k 2 I°^)ldtf>Bmh(T(z), (2.64) 

d(j) ^ cosh a(z) — \/o(z)dt sinh a(z). (2.65) 

This boost representation is known for the Myers-Perry black hole [3, 6]. However our lead¬ 
ing order solution covers much wider class of solutions including the asymptotically AdS 
black holes. This boost representation is the universal feature for the large D stationary 
black holes under our ansatz. 

The boost parameter, cr(z), is not constant and the boost transformation is inhomoge¬ 
neous in z-direction. This z-dependence also appears as the horizon radius inhomogeneity 
in g pp . Thus the identification with the Schwarzschild black hole by the boost transforma¬ 
tion is valid only locally in ^-direction. However this inhomogeneity along z-direction is 
not so crucial since the dynamics along ^-direction is sub-dominant compared with the ra¬ 
dial dynamics. This boost transformation is the property of the radial dynamics, which is 
the leading order dynamics of large D black holes. The local effect by z-dependence would 
be introduced at the sub-leading order in 1/D expansion as corrections. Such z-dependent 
dynamics in the higher order structure of the 1/D expansion contains the essential effects 
by the rotation and horizon topology. 

3 Perturbation 

As one interesting application of the boost property in the leading order solution we 
consider the perturbation problem. By using the boost transformation we can obtain 
the solution of the perturbation equation of the stationary black hole from one of the 
Schwarzschild black hole. In the following we investigate the decoupled mode perturbation, 
whose the frequency is u = O(D 0 ) [7], of the asymptotically flat black hole given by 

V 0 (z) = 1, (3.1) 

in the leading order metric in section 2. The extension to the asymptotically AdS black 
hole is straightforward. 
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3.1 Schwarzschild black hole 


We study the perturbation on the Schwarzschild black hole. Next we perform the boost 
transformation on the perturbation solutions to obtain the perturbation solution of the 
stationary solution. There are various ways in the representation of the perturbation 
around the Schwarzschild black hole [27, 28]. Here we fix the gauge for the perturbation, 
hp U , to solve the perturbation equation. Especially we use the transverse traceless gauge 
as 


V = 0, gT V = 0, (3.2) 

where is the covariant derivative of the D dimensional background metric g^ u . Then 
the perturbation equation becomes equivalent to the Lichnerowicz equation given by 


V p v p V + 2 vr = o, (3.3) 

where Rp U pa is a background Riemann tensor. The explicit components of the perturbation 
equation are complicate so we do not show them here. The background spacetime is the 
D = n + 3 dimensional Schwarzschild black hole given by 


ds 2 = —f(r) 1 dt 2 + f(r) 1 dr 2 + r 2 (d6 2 + sin 2 Odcf 2 + cos 2 9dQ%_ 4 ) , (3.4) 


where 


fir) = 1 



(3.5) 


Here we introduced the new radial coordinate R = ( r/rS ) n . This radial coordinate R is 
related with the radial coordinate p which we used in the previous section by 


R = cosh 2 p. 


(3.6) 


The perturbation is decomposed into the scalar, vector and tensor type with respect to 
S D ~ i in the metric 5 . In this paper we investigate the scalar and vector type perturbation 
on S D ~ 4 . Using the vielbein 

4°ch = 4 W)dt , egch = —r ==w> 4ch = rde ’ 41 =rsin9d(t>, = r cos<?e (A) , (3.7) 

v/( r ) 

where e^ A ' ) is a vielbein on the unit sphere S D ~ 4 , the scalar type perturbation, h\fj, on 
can be written by 


h( s )(j T rd T u — [ f( s )S cll Y( s ) p (h SR RS'lSdiQ y(S) (J) (A) 

n pv ax ax - e e J/j e sch e Sch + Ji e Sch e Sch 


+ r 2 cos 2 e(H^ S)Sch Yf- ] 1AB + H { L S)Sch Y ( A ) B j )dx A dx B ], (3.8) 


5 Note that the usual scalar, vector and tensor type perturbation of the Schwarzschild black hole is 
based on the decomposition on S D ~ 2 [27]. Then, for example, our scalar type perturbation on S D ~ 4 
consists of the scalar, vector and tensor type perturbation on S D ~ 2 . 
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(3.9) 


where /, J = 0,1, 2, 3. Y ^ is the scalar harmonics on S D 4 satisfying 

(D 2 + Xs)Yf = 0, 

where the eigenvalue Xs is given by 

As =j(j + D- 5). 

(S)i 

The scalar derived tensor harmonics Y AB is defined by 

yJS 7 ' = D^DsYf } + -^Y f>q AB . 

The vector type perturbation, h^\ on S D ~ 4 is given by 

fj V)Sch Y { A )j e^e^ A) + R (V) Y ^ e (^) e (s) 


h^dx^dx" = 


where Y^ 2 and Y^fn are the vector harmonics on S D 4 defined by 


l A ai m ii AB 

>.2 


(D 2 + Ay)Y^' = 0, D a Y^"=0, Y^ = D (a Y b) 
with the eigenvalue Ay given by 

Ay = j(j + D - 5) - 1. 


^A w (Y)j 


AV)j 


(Y)j 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


The perturbation equation is a coupled PDE system for scalar and vector type pertur¬ 
bations. At large D the situation is changed and the perturbation equations become 
decoupled ODE system. We can find decoupled perturbation variables p^ Sch for the 
scalar type perturbation easily as 


F 


(S)Sch _ „(S)Sch 


_ . r 

00 ~ J 00 + J 11 


(S)Sch p (S)Sch _ AS )Sch „(S)Sch p (S)Sch _ AS )Sch 

j -foi — Zoo Zn i -^oi — Jo 


rp{S) Sch _ AS) Sch p (S)Sch _ ,(S)Sch p (S)Sch _ AS)Sch p (S)Sch _ „(S)Sch 

^02 ~ J 02 i ^03 — Z 03 ) ^ 12 — Z 12 J ^13 — Zl 


/01 

f (S 

M3 


K 


(S)Sch _ „(S)Sch ,-,(S)Sch _ „(S)Sch F (S)Sch _ AS )Sch 


22 


= /: 


22 


) ^ 2.3 


= /: 


23 


77’V 

5 -^33 


= /; 


33 


(3.15) 


and 


F 


(S)Sch 


= /} 


(S)Sch 


ri(S)Sch 

T,L 


= Hr 


(S)Sch 


T,L 


(3.16) 


The vector type perturbation has following decoupled perturbation variables at large D 


(y)Sch _ f (v)Sch F (y)Sch _ ^-(v)Sch 


= /, 


(3.17) 


Then each decoupled perturbation variables are expanded at large D as 


F, 


(S)Sch 


(k) p(S)Sch 


(k) j )S c h 


ab 


E ± ab _ p(S)Sch _ \ ' ra _ p,(S)Sch _ \ ' 

B k ’ a 2 —j r,k ’ T,L 


fc >0 


k> 0 




(k) p(S)Sch 

T ’ L (3.18) 


k> 0 


n * 
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for the scalar type perturbation and 


F (V)Sch = y- (4 j?j y ) Sch ^ p (V)Sch 
k> 0 


E 


(k) p(y)sch 


n*■ 


(3.19) 


for the vector type perturbation. By using transverse traceless gauge conditions we can 
obtain F^ Sch . F^ Sch and ' |Sdl from F^'p Scl ' and f\ X iSch . Thus the equations to be 


solved are the equations of Fjj’ Sch for the scalar type perturbation and F^ ^ Sch for the 
vector type perturbation. The perturbation equations for are 


^ R(R 


d_ 

dR 

d_ 

SR 

d_ 

dR 


R(R 
R(R — 
R(R — 


a m F (s >sch w 4 s)Sch 

"sR ' n R -1 

_ ( * )J o? )Ml 

’dR 01 R — 1 


4 R (R 

dR v 


d_ 

dR 


R(R — 


Ji R < R 


_ (t) 4 S)Sch 

’dR 02 4R(R — 1) 

Amw* _ w 4 s)sdl 

j <9R 03 4R(R — 1) 

.. a m (s)sch WF<f )Sch 

’flR 11 R 

,,«(»„» _ w 4 )Sdl 

j <9R 12 4R(R — 1) 

Aysia _ (t) 4 s)sdl 

j <9R 13 4R(R — 1) 

— RfR — i)F’ S ’ Sch 
5R K K 1 <9R ^ 22 

— RfR — i)—( k '>F ( ' S ' >Sch 
3R 1 ’dR 23 


(k) AS) 
°00 ) 

(fc) 4?, 


(fc) c(S) 
°02 > 


(k) AS) 
°03 > 

(fc) c(S) 

°13 > 

<‘>51?, 

: <‘>S<f. 


(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 


The k —th order source term is coming from the lower order solutions and the 

source term does not have the decoupling property. The perturbation variables, ( k )Fjj’ Sch , 
depend on R and 6 6 . However the perturbation equation becomes ordinarily differential 
equation with respect to R at large D. Furthermore each fc-th order variables ( k )Fjj^ ch 
are decoupling and we can solve the solution by the straightforward integration procedure 
at each order. The vector type perturbation has the same property. The perturbation 

6 The Schwarzschild black hole has a spherical symmetry. So the perturbation equation can be reduced 
to the ordinarily differential equation in principle. Here, for the usefulness in the following discussion, 
we use the harmonics on S D ~ 4 and the perturbation equation becomes the apparent partial differential 
equation. 
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equation for ' )Sch is 


9 P(P 1) 9 ( fc )F (y)Sch {k)F " )&dl - ^S {V) 

8R [ j 9R 0 4R(R — 1) _ ’ 

(3.30) 

%<r i) 9 (t>F ! V>Sch _ m s m 

9R 1 j 9R 1 4R(R —1)“ 61 ’ 

(3.31) 

AR(R-i)F(^y>^ = m s p, 

(3.32) 


(3.33) 


The leading order solution satisfying the regularity condition on the horizon and at infinity 
is given by 


(0) P (S) Sch _ A(0) (0) (S)Sch _ A(0) (0) (S)Sch _ B(0) 

00 R — 1 ’ 01 R — 1 ’ 02 V /R(R_ i )’ 

(0) p(5)Sch _ c(6) ( 0 ) (S) Sch _ D{0) (0) p(S)Sch _ B{6) 

03 V /R(R- 1)’ 11 R ’ 12 V /R(R — 1) ’ 

(0) i7i(>S)Sch _ C'(O) j7>(S)Sch _ n (o) jpfSISch _ n (o) ^(SlSch _ n 

-*13 — V /R(R _ xj ’ ^ 22 _ U ’ ^ 23 _ U ’ ^ 33 _ ' 


(3.34) 

(3.35) 

(3.36) 


for the scalar type perturbation and 


(Q^tOSch 


W ( 6 ) fOlriWSch W(0) fo) ^(y)Sch ( 0 )^(V)Sch n fo 

^R(R-l)’ 1 “VRCR-i)’ 2 “ 3 ,( ’ 


for the vector type perturbation. The leading order perturbation solution contains un¬ 
determined integration functions A(0), B(0), C(0), D(0) and W(0) for each type pertur¬ 
bation. If we go to the higher order of k we can obtain the non-trivial relation between 
integration functions. Especially the solution at k = 1 and k = 2 order gives quasinormal 
mode frequency of the Schwarzschild black hole for the vector and scalar type perturbation 
respectively. The detail analysis and results of this will be given in section 4 and Appendix 
D. 


3.2 Stationary solution 

We apply the boost transformation to the perturbation solution of the Schwarzschild 
black hole. At first, before taking large D limit, we give the general procedure for the 
perturbation of the stationary black hole solution given by 

ds 2 = —A{p, z) 2 dt 2 + B(p.z) 2 dp 2 + F(p, z) 2 (dc/) — w(p, z)dt) 2 

+ G(p,z) 2 dz 2 + H{p,z) 2 dVL 2 D _ A . (3.38) 

We consider the scalar and vector type perturbation on S D ~ 4 . To give perturbation 
variables we use the following vielbein 

e^=A(p,z)dt, = B(p,z)dp, = F(p, z)(d<j> — w(p, z)dt) (3.39) 
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and 


e (3) = G(p, z)dz, e (A) = H(p,zp A \ (3.40) 

where eh) is the vielbein on the unit sphere S D ~ 4 again. Then the perturbed metric h /J:iy 
is given by the components with respect to these vielbeins as 

hffldx^dx 1 ' = e ~ iujt Yfe^e^ + /{ S) D A ¥f ) e (/) e (A) 

+ H(p, zfiH^YfqAB + H ( f ) Y^ ) J)dx A dx B '\ , (3.41) 

for scalar type perturbation and 

h ( y)dx> x dx v = + H { pY^e^e^ , (3.42) 

for the vector type perturbation. Then we perform the boost transformation on these 
perturbed metrics of the leading order solution obtained in section 2. The boost transfor¬ 
mation gives the relation between vielbeins as 

2 — ro(z) 2 (l — sinh 2 p) 


D (°) 

"Sch 


= (0) 


— sinh p 


D ( 2 ) 

-Sch 


ro(z) 2 coshpy 1 + r 0 (z) 2 sinh 2 p 

coshp ( 2 ) 2 + ro(z) 2 sinh 2 p 

^1 + ro (z) 2 sinh 2 p r ° ^ 


1 - r 0 (z) 2 
1 + r 0 (z) 2 sinh 2 p 



d 2 ) 


1 - r 0 (z) 2 


M 


r 0 (z) sinh 2 p 


(3.43) 


(3.44) 


where we used the vielbein of the Schwarzschild black hole (3.7) and of the stationary 
solution (3.39) for the leading order metric (2.61) ' . This relation between vielbeins gives 
also the relation between the perturbation variables of the Schwarzschild black hole and 
general stationary solutions. Since we know the decoupled perturbation variables on the 
Schwarzschild black hole, we can obtain the following decoupled perturbation variables, 
Fj’V , Fj S \ Fpj. f\ X ^ and F^ ^ of the general stationary solution at large D as 


*00 — 


cosh 2 p 


1 + ro(z) 2 sinh 2 p 


AS) 2 sinh P cosh Py/ 1 ~ r o( z ) 2 AS) 

J00 1 ■ ,2 J 02 


1 + ro(z) 2 sinh 2 p 


+ 


(1 ^ r 0 (z) 2 )sinh 2 p f{S ) 


1 + ro{z) 2 sinh 2 p 


fh’ + f 


C s) 

n > 


(3.45) 


p(S) _ 
-*11 — 


cosh 2 p 


1 + ro(z) 2 sinh 2 p 


AS) _ 2 sinhpcoshpy 7 ! - r 0 (z) 2 (s) 
00 1 + ro(z) 2 sinh 2 p 02 

(1 -r 0 (z) 2 )sinh 2 p {S ) _ (S) 

1 + ro(z) 2 sinh 2 p 22 11 


(3.46) 


^(S) _ cosh 2 p - (1 - r 0 (z) 2 ) sinh 2 p f ( S ) 
-CQ2 — • i 2 J 02 


1 + ro(z) 2 sinh 2 p 

\/l - r oi z ) 2 sinh p cosh p^ as) , 

1 , M2 • u2 WOO + 9 02 h 

1 + ro(zy smh p 


We consider the asymptotically flat black hole Vo (z) = 1. 


(3.47) 
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and 


p(S) _ 

r 22 ~ 


cosh 2 p 


AS) 2sinhpcoshp v / l - r 0 (» 2 f ( S ) 

2 J22 1 , „:„i,2 „ j 02 


1 + ro(z) 2 sinh 2 p J22 1 + ro(z) 2 sinh 2 p 

(1 -r 0 O) 2 )sinh 2 p f[s) 

Jo 


+ 


1 + ro(z) 2 sinh 2 p 


2 J oo > 


(3.48) 


p(s) _ 

r 01 — 


cosl1 '’ ^f-sinhpj , , 1 -;° ( ; ) 1 2 /<f, 


1 + r 0 (z) 2 sinh 2 p 


1 + vq{z) 2 sinh 2 p“ 


(3.49) 


p(s) _ 

r 12 ~ 


C0ShP /Sf>. (3.50) 


1 + tq(z) 2 sinh 2 p 


1 + r 0 ( 2 :) 2 sinh 2 p“ 


p( s ) _ 

^03 ~ 


cosh p 


1 + ro(z) 2 sinh 2 p 


=/o ( 3 } ~ sinh p a 


1 - r 0 ( 2 :) 2 


1 + r 0 (z) 2 sinh 2 p 


f: 


(S) 


2 2 23 > 


(3.51) 


^23 — 


cosh p 


1 + ro(z) 2 sinh 2 p 


=/23 } ~ sinh p * 


l~r 0 (2:) 2 j( g) 


1 + ro(z) 2 sinh 2 p 


2 /o3 > 


(3.52) 


Ff = /r-sinhp. - jf>, 


1 + ro(z) 2 sinh 2 p 


1 + ro( 2 :) 2 sinh 2 p “ 


(3.53) 


_p( s ) — 
r 2 — 


cosh p (5) / 1 - r-o(^) 2 ,(S) 

= / 2 sinh/^/ 1 f __ , 2 .,2 4 0 > 


1 + ro(z) 2 sinh 2 p 


1 + ro(z) 2 sinh 2 p“ 


(3.54) 


p(S) _ f(S) p(S) _ f (S) P (S) - f( 5 ) p( 5 ) - f( S ) p( s ) - w( 5 ) 

-*13 — 213 i ^33 — 233 > r l ~ Jl i r 3 ~ 23 j r T,L ~ 


(3.55) 


for the scalar type perturbation and 


- 


coshp AV)_. h / 1 - rp{z) 2 (V) 

=/ 0 smhpy i | __ t 2 Ji » 


y/1 + ro( 2 () 2 sinh 2 p 


1 + vq(z) 2 sinh 2 p 


(3.56) 


F (V0 _ 
r 2 — 


cosh p 


yj 1 + ro(z) 2 sinh 2 p 


/ 2 (V) -sinhp. 


f ~ l’o(^) 2 ,(V) 


1 + ro(z) 2 sinh' p 


2 20 5 


(3.57) 


and 


(V) _ AV) AV) _ ,00 


TJi\ v J f\ v J TP\ V ) _ f 

*1 — 21 ) ^3 — 2; 


F (^) _ ^(V) 
3 1 r L ~ L ) 


(3.58) 
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for the vector type perturbation. Each perturbation variables are expanded at large D as 


and 


r (s) r (s) r is) 


fc >0 


fc >0 


n" 


T.L 


fc >0 


n" 



MfW W fP 

/ -j ^k 5 T,L / j^k 

k> 0 /c>0 


(3.59) 


(3.60) 


Then the perturbation equation for these decoupling variables can be obtained from the 
Lichnerowicz equation under the transverse traceless gauge condition. Now we have ob¬ 
tained the leading order metric of the stationary black hole in the l/D expansion in section 
2. Then the perturbation equation for k = 0 order is also obtained by using the leading or¬ 
der metric. The perturbation equation of the decoupling variables can be found to satisfy 
same equations with one of the Schwarzschild black hole at the leading order by changing 
the radial coordinate from p to R by 8 


R = cosh 2 p. 


(3.61) 


At higher orders the perturbation equation of the stationary solution becomes different 
from the boosted one of the Schwarzschild black hole in source terms. This is because 
the boost relation is valid only in the leading order metric, and higher order structure 
gives the essential effect coming from the difference from the Schwarzschild black holes. 
However the homogeneous part of the equation at each order is always same with one of 
the Schwarzschild black hole. So, in principle, the perturbation equation is a set of decou¬ 
pled ordinarily differential equations and we can solve it by a straightforward integration 
method with the Green’s function. 

The regular leading order solution of the perturbation equation for the stationary black 
hole is 


K — 1 K — 1 yj 

C ( 9 ) (o) jp(S) _ D{0) (0) p(S) = _ 

R(R — 1) ’ 11 R ’ 12 V 


for the 


(o) F (S) = C{6) 

03 a/R(R — 1) ’ 

(0) 4 S) = (0) 4 S) = o, (0) 4 S) = 0, (( 

scalar type perturbation and 

(0) F (V) = W{6) (0) F (V) = W(0) 0) F (V) = (0) F 

0 R i R 


(3.62) 

(3.63) 

(3.64) 


(3.65) 


for the vector type perturbation. In the derivation of the perturbation equations and 
solutions we need only the leading order solution of general stationary solutions given in 

8 This relation can be understood by the horizon position. The Schwarzschild black hole has the horizon 

at R = 1 and the general stationary black hole horizon position is set to be p = 0. 
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section 2. To go to more higher order perturbation equations, the higher order solutions 
of stationary solutions are required. In the next section we consider the higher order 
perturbation equation of the Myers-Perry black hole, which is a known exact solution. 

Note that the quasinormal mode frequency of the decoupled mode perturbation is ob¬ 
tained by solving higher order perturbation equation in 1/D. The higher order structure 
of the perturbation equation reflects the higher order structure of the background metric, 
which does not have the boost property. Hence the frequency does not have the boost 
property either. This implies that the quasinormal mode frequency of the Myers-Perry 
black hole cannot be obtained by the boost transformation from one of the Schwarzschild 
black hole. This non-boost property allows the existence of the instability mode in the 
decoupled sector of the perturbation of the Myers-Perry black hole. In contrast the quasi¬ 
normal mode of the non-decoupled sector perturbation is obtained in the leading order 
structure and we expect that such mode has the boost property as the universal feature 
[5], 


4 QNMs of Myers-Perry black hole 


We pursue the decoupling property of perturbation variables at large D expansion shown 
in the previous section to more higher order structure to obtain the quasinormal mode 
(QNM) frequency of the D = n + 3 dimensional singly rotating Myers-Perry black hole. 
This was calculated numerically in [29] for non-axisymmetric perturbations. Here we 
give the analytic formula of the QNM frequency both for the axisymmetric and non- 
axisymmetric perturbations. The metric of the Myers-Perry black hole is given in the 
Appendix A and it is written by the following parameter setting 

V 0 (z) = 1, «=^, n H = 2 > = cos0 ' r o(0)~ 2 = 1 1 X a 2a , (4.1) 

2 1 + a z 1 + a z cos z u 

in the leading order metric obtained in section 2. Our detail analysis in this section con¬ 
centrates on the scalar type perturbation since the analysis of the vector type perturbation 
is performed in the same manner. The only result of the QNM frequency for the vector 
type perturbation will be given later. The scalar type perturbation variables are defined 
in the same way with one in section 3 as 


h^dx^dx" = e - iujt e im * 


LbYfe^e^ + f a D A Y) J> e^e' 


7 ^)Ja)(A) 


+ r 2 cos 2 e{H T Y { pq AB + H L Y ( f ) B j )dx A dx B ], (4.2) 


where the vielbein was defined in section 3 9 . So the perturbation is parameterized by 
w, m and j. j describes the deformation parameter on S D ~ A . Thus j = 0 correspond 
to the perturbation without the deformation on S D ~^. Then the decoupled perturbation 

9 We omit the index ^ representing the scalar type perturbation in section 3. 
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variables can be constructed and expanded at large D as 


Fu = 


f: 


(k) 


u 


Fa = ^~^' = 


(*) 


k> 0 


n' 1 


fc >0 


k> 0 


p( fc ) 

r T,L 

U ' 

n K 


(4.3) 


where D = n + 3. The boundary condition for the perturbation on the horizon is the 
ingoing boundary condition [29] as 


0) = KAr, 0)(r - 2k , 


(4.4) 


where r+ is the horizon position and k is a surface gravity. The component of h^ u in the 
Eddington-Finkelstein ingoing coordinate should be regular on the horizon. Here we omit 
the t, cj) and x A dependence just for the simplicity 10 . The boundary condition at infinity 
R 1 is the decoupled mode condition [7] as 


V = 0 ( R- 1 ). 


(4.5) 


This condition is equivalent to the outgoing wave boundary condition at infinity for the 
decoupled mode defined by uro = 0(D°). 

In the following we give the brief summary of perturbation solutions of the scalar type 
perturbation up to k = 2 order. At k = 2 we obtain QNM frequency for the first time. 


Leading order solution (k = 0) The leading order solution satisfying the boundary 
conditions is 


p( 0) A (0)( 6> ) F (0) ^4(o) (^) 

00 R - 1 ’ 01 R - 1 ’ 



C(Q)(0) p(o) 
v/RO^T)’ 11 


F (o) C(o)W 

13 v/R(R-l)’ 


D(o)( 9 ) 

R 


51 ( 0 ) 

22 


= 0, 




v/R(R-i)’ 
B { 0 )(g) 
v/R(R-i)’ 


f £ ] = 0, 4 0) = 0. 


(4.6) 


At leading order we have four independent integration functions, which cannot be deter¬ 
mined by the boundary condition. The solution for Fj 0 ^ and F^ L can be obtained by the 
transverse traceless gauge condition, and we omit them since it is not so important. 

10 x A is a coordinate on S D ~ 4 . 
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Next to leading order solution ( k = 1) The next to leading order solution of Fjj is 
represented by 


771(1) ^(1) ( d ) p( l) 

^00 - “—— + tv 


R- 1 


00 1 01 


F <:> = Mfl+ E<i> 


p(i) _ 

^02 — 


Bg )W 

v/R(R-i) 


_i_ aU) p(i)_ 

+ -^02 5 t ' Q3 — 


17 ( 1 ) _ jD (l)( 0 ) , £,( 1 ) 
^11 - —--l-F 


R 


pW _ 

11 > r 12 — 


R- 1 

C { i)(0) 
\/R( R ~ 1 ) 
fl(p(g) 

v/R(R -1) 


01 j 


, p(i) 

+ -^03 ) 


+ 

+ r 12 > 


p(i) _ 

-^13 ~ 


c w (0) 
v/R(R-i) 


+ F^ F^ — F^ F^ — 0 F^ — F^ 

^ r 13 5 r 22 — r 22 > ^23 — u > ^33 — r 33 > 


(4.7) 


where the solution with the hat comes from the integration of the source term composed 
by the leading order solution. At this order we have new four integration functions Ar^{6), 
C(ij(6) and D^(0). To satisfy boundary conditions we obtain some non-trivial 
relations for A( O )(0), F( 0 )(0), (7(0) (0) and D^(9). Conditions we find are 


B {0) (9) = 0, (7 ( o)(0) = 0, 


(4.8) 


and 

cos 0(1 + a 2 cos 2 0)^(1 + a 2 ) cos 0A( O ) (0) — sin OD 1 ^ (0)^ 

+ [j + (—1 + iam + a 2 (—3 + 2 j — iui) — iui) cos 2 0 

+ a 2 (l + iam + a 2 (j — 1 — iui) — iu) cos 4 0] F>( 0 ) (0) = 0. (4.9) 

Using this condition we can eliminate A( O )(0). To obtain QNM frequency we should specify 
the function F>( O )(0), and we can do it at next to next to leading order. 

Next-to-next-to leading order solution (k = 2) At this order we can obtain the 
solutions in the same way with the next to leading order. We find additional non-trivial 
condition for integration functions. By eliminating A( O )(0), the additional conditions are 
written as 


Cl (0)^") (0) + c 2 (0)C ( ' 1} (0) + c 3 (0)C (1 ) (0) 

+ c 4 (0)F[ / o) (0) + c 5 (0)F; o) (0) + c 6 (0)F (o) (0) = 0, (4.10) 


and 


d!(0)Ff' } (0) + d 2 (0)Ff o) (0) + d 3 (9)D{ 0) (9) 

+ d 4 (0)F( O )(0) + d 5 (0)C { , 1) (0) + d 6 (0)C' ( i ) (0) = 0. (4.11) 

Each functions Cj(0) and dj(0) have messy forms so we do not show them here. They 
are given in the Appendix E. In general it seems that we should solve these differential 
equations to obtain the QNM frequency. However, we find that the QNM frequency can 
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be obtained only from the boundary conditions 11 . At 9 = 7t/ 2 the boundary condition 
is read from eqs (4.10) and (4.11) as 


C'(i)( 6 ')l( 9 = 7 r /2 = cos - 7 3 6 >, D( O) (0) | e=7r/2 = cos - 7 0. (4.12) 

On the other hand the behavior at 9 = 0 is 

C {1) (9)\ e=0 = sin* 5 " -1 9 , D {0) (9)\ e=0 = sin 5 " 9 , (4.13) 

where 6 ^ depends on a, m, j and uj. To see which value 5^ should take, we observe the 
spherical harmonics, E> e (9)e im ^Yj S \ on S D ~ 2 , defined by 

[A s d-2 + 1(1 +D- 3)] S*e im *Y^ s) 


■ 9 1 D-Aft lft ( sin6>cosD 4 4 S< 
sin 0 cos u 4 # dO \ a9 


” 2 ;S ' - s' + ^ + D - 3)S< 


e im0 Y (S) = Q 


sin" 9 cos 2 9 

The solution of this equation can be written by the hyper geometric function as 
= cos- 7 9 sii:J m l 9 2 F 1 (—ks, ks + \m\ + j + n/2, \m\ + 1; sin 2 9), 


(4.14) 


(4.15) 


where ks is a non-negative integer satisfying 

l = j + \m\ + 2ks- 


(4.16) 


At large D this solution becomes 

= sin l - j 9 cos j 9 + 0(1/D) 

= sinl m l +2fcs 9 cos- 7 9 + 0(1/D). (4.17) 

ks describes the overtone number of along 9 direction, and this overtone number can 
be observed from the behavior of at 9 = 0 at large D limit as seen eq. (4.17). From 
this observation on the spherical harmonics we impose the following condition on <5 ^ as 
the harmonics condition 


&t=l~3 

= 2 k s + \m\, (4-18) 

where l is parametrized as eq. (4.16) by non-negative integer ks■ Then eqs. (4.10) and 
(4.11) under the condition (4.18) give one non-trivial algebraic condition on uj. This non¬ 
trivial condition on uj can be regarded as the QNM condition. In the following we solve 

11 Actually we can confirm that eqs. (4.10) and (4.11) can be solved regularly under the QNM frequency 
obtained below. Thus there may be some non-trivial structure in the solutions of eqs. (4.10) and (4.11) 
and the boundary conditions. 
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the condition eq. (4.18) for oj in some cases. The explicit form of the QNM condition 
(4.18) is given in the Appendix F. 

Note that the condition (4.18) is equivalent to the regularity condition to derive eq. 
(2.56) in section 2. The condition (4.18) is suggested by the analysis for spherical har¬ 
monics at large D. However, as we can see in Appendix C, the analysis of the spheroidal 
harmonics also gives same condition (4.18). Thus we expect that the condition (4.18) 
would hold also in the gravitational perturbation of the Myers-Perry black hole. 

4.1 QNM frequency 

We show some explicit results of the QNM frequency of the scalar and vector type pertur¬ 
bations by solving the QNM condition (4.18). The instability exists only in the scalar type 
perturbation and the vector type perturbation is always stable for the decoupled sector 
perturbation 12 . We set the horizon radius to unity as 

r 0 = 1, (4.19) 


by fixing the unit. 

4.1.1 Scalar type perturbation 

The QNM frequency is given by the solution of the algebraic equation (4.18), which can be 
solved analytically. We solve eq. (4.18) for some cases below. Since the j / 0 modes do not 
show any instability, we consider only j = 0 modes. The perturbation with j ^ 0 describes 
the deformation of S D ~ 4 of the Myers-Perry black hole. So we might be able to say that 
the axisymmetric instability occurs only in the ”S-wave” sector in the perturbation similar 
to the Gregory-Laflamme instability of the black brane [32, 33]. 

Then, for this mode, the quasinormal mode condition (4.18) can be written in a rela¬ 
tively simple form as 

(1 + a 2 ) 3 w 3 — (1 + a 2 ) 2 (4i — 3£ + 3 am)u 2 

+ (1 + a 2 ) (—4 + 7£ — 3£ 2 — 6 iam(£ — 1) + a 2 — 4 + 3£ + 3 m 2 )oj 
+ (—i£ 3 + £ 2 (3 i + ia 2 + 3 am) — £(2i + 5 am + a 3 m — ia 2 (—2 + 3m 2 )) 

+ am( 2 — 2 iam — a 2 (m 2 — 2))) = 0. (4.20) 

At a = 0 the scalar type perturbation has the physical degree of freedom for £ > 2 [30] 
and, hence, we study the modes of £ > 2. 

12 The instability mode is conjectured to exist only in the decoupled sector (saturated mode) [29]. The 
tensor type perturbation of the singly rotating Myers-Perry black hole does not have the decoupled sector 
and it was shown to be stable [31]. 
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Figure 1: The QNM frequency for ( £,m,j ) = (4,0,0). The black thick line is the pure 
imaginary mode and it shows the instability at a 2 > 3. The gray thick line is the stable 
mode. 


Schwarzschild black hole At a = 0 where the solution is the Schwarzschild black hole 
we can solve eq. (4.20) explicitly by 

uj = ±y/£-l-i{l-l), oj = -i£. (4.21) 

These modes correspond to the decoupled scalar and vector type mode of Schwarzschild 
black hole on S D ~ 2 [7]. 

Axisymmetric perturbation Next we consider the axisymmetric perturbation m = 0 
with a/0. In this perturbation we find that there are stationary perturbations uj = 0 at 
a c where 


a 2 c = £-1 

= 2 k s - 1, (4.22) 

and we used eq. (4.16) with j = m = 0. The assumption of £ > 2 implies k$ > 1 
for j = m = 0, and this result reproduces the result (2.56) derived from the effective 
theory obtained in section 2. In figure 1 we give the plot of the QNM frequency with 
(. £,m,j ) = (4,0,0). The black thick line, purely imaginary frequency mode, shows the 
instability mode and it becomes the marginally stable mode at a 2 = 3. For the axisym¬ 
metric perturbation with uj = 0 and j = 0 we can solve eqs. (4.10) and (4.11) explicitly 
by 

~ sin 1+a 0 

Cb(0) = O, D o (0) = D o — -2-(4.23) 

1 + cos z 0 

where Dq is a constant. Then we can see that this stationary solution is equivalent to the 
perturbation solution (2.55) obtained in the effective theory. This equivalence suggests 
the validity of our QNM condition (4.18). 
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Figure 2: The plot of the bar mode QNM frequency of ( £,m,j ) = (2,2,0) is shown. The 
black dashed line is the superradiance condition. The onset of the dynamically unstable 
mode is same as the one of the superradiance instability. 


Non-axisymmetric perturbation The non-axisymmetric mode perturbation also shows 
the instability at a 2 > a 2 where a c is given in eq. (4.22) as 

a 2 c = £-1 

= \m\ + 2k s — 1. (4.24) 


At the critical rotation the instability mode satisfies the superradiance condition as 

a c m 


OJ I n —a,„ - 


1 + at 


= m£l H {ac)- 


(4.25) 


Hence the dynamical instability mode shows the superradiance instability at the same 
time. This coincidence of the onsets of dynamical and superradiant instability has been 
also observed in the QNMs of the Myers-Perry black hole with equal spin [6] up to 0(1/D) 
correction. The numerical results [34] also suggest the coincidence for the equal spin case. 
On the other hand the numerical analysis of QNMs of the singly rotating Myers-Perry 
black hole does not show such coincidence [29, 35]: the superradiant instability appears at 
slower rotation than the critical rotation of the dynamical instability. Thus our coincidence 
is just the property only at large D limit. If one considers 1/D corrections of QNMs of 
singly rotating Myers-Perry black hole, there would be a difference between the onset of 
the superradiant and dynamical instability. 

For the bar mode defined by m = £, we can solve eq. (4.20) explicitly by 


±\/m — 1 + a(m — 1) m — 1 T ay/m — 1 
oj± =-——^- i- 


1 + a 2 * 1 + a 2 ’ W(0) 1 + a 2 *1 + a 2 ' 

cj + shows the instability at a > a c and saturates the superradiance condition at a = a c . In 
figure 2 we show the plot of the bar mode QNM frequency for m = 2. The QNM frequency 
for the bar mode was obtained numerically in [29]. The behavior of the frequency shows 
good agreements with our analytic results. 

The non-axisymmetric perturbation with l > rn also shows the instability. But its 
origin is different from one of the bar mode. The instability mode of the bar mode 


am 


m 


— i- 


(4.26) 


27 


















Re[co»r+] 



Figure 3: The plot of the QNM frequency 
dashed line is the superradiance condition, 
while the black thick line shows the unstable 
exactly when the superradiance condition is 
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of ( £,m,j) = (4,2,0) is shown. The black 
The gray thick line is the stable mode w_, 
mode w-|_. The dynamical instability appears 
satisfied. 


perturbation comes from the scalar type perturbation of the Schwarzschild black hole 
at a = 0 where the frequency is complex (see eq. (4.21)). On the other hand non- 
axisymmetric mode with £ > m has an instability in the mode coming from the vector 
type perturbation of the Schwarzschild black hole, whose frequency is a pure imaginary as 
(4.21). This is same situation with the perturbation of the Myers-Perry black hole with 
equal spin [6]. In figure 3 we give the plot of QNM frequency for £ = 4, m = 2. As we can 
see, the mode showing instability becomes purely imaginary at a = 0, which corresponds 
to the vector type perturbation of the Schwarzschild black hole. 


4.1.2 Vector type perturbation 


The QNM frequency of the vector type perturbation can be given as 


^ = a; V 0 + ^NLO +0(1/rt 2 ) 
n 


where 


and 




am 
1 + a 2 



t~3 \ 
1 + a 2 / 


w nlo = ^^2 (^ - 1 + ° 2 0' “ !) + log (1 + a 2 )) 

_ i f(j _ 1)2 + (! + a 2 )( l - j)(e + 3-2)- 2 a 2 m 2 


(1 + a 2 ) 2 

. . , l-j 2(1 — j) 


1 + a 2 (1 + a 2 ) 2 


log (1 + a 2 ) . 


(4.27) 


(4.28) 


(4.29) 


We have solved the perturbation equation up to k = 2 order and obtained the QNM 
frequency up to 1 /D correction. This QNM frequency is obtained by imposing the same 
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harmonics condition with the scalar type perturbation on the behavior at 6 = 0 motivated 


by eq. (4.17). The QNM frequency (4.28) reduces to the vector type perturbation on 
S D ~ 2 of the Schwarzschild black hole [7] at a = 0 as 



(4.30) 


5 Summary 

The large D expansion method has been found to be useful to solve the gravitational 
problem as shown in [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. In this paper we have performed the 
next step for the further developing of the large D expansion method by constructing 
the effective theory of the large D stationary black hole in asymptotically flat or AdS 
spacetime. Considering the additional degree of freedom of the horizon, rotation radius of 
the black hole, the effective theory becomes non-trivial and allows the existence of various 
solutions compared to the static case in the asymptotically flat spacetime which has only 
one unique solution, Schwarzschild black hole [9]. This new degree of freedom appears as 
the Lorentz boost of the black hole along the rotational direction, and we can rewrite the 
solution of the effective equation in very simple and illuminative expression by using this 
boost property. 

As applications of the effective theory we have considered the ellipsoidal embeddings 
and obtained the Myers-Perry black hole and bumpy black hole solutions in perturbative 
manner. Then we have succeeded to derive the threshold angular momentum of the insta¬ 
bility of the Myers-Perry black hole. The bumpy black hole construction and perturbation 
analysis of the Myers-Perry black hole had been investigated numerically. These numerical 
studies of the solution and perturbations need very advanced and sophisticated technique. 
However the effective theory obtained in this paper is rather simple equation and we can 
obtain the solution easily. 

As another application we have obtained the quasinormal mode frequencies of the 
singly rotating Myers-Perry black hole. Our effective theory is only for the stationary 
black hole. But, if we consider the dynamics with wro = O(D 0 ), the time-dependence 
becomes sub-dominant compared with the radial dynamics. As a result we could solve 
the perturbation equation as the ordinarily differential equation with respect to the ra¬ 
dial direction. Furthermore the leading order solution of our effective theory has a very 
important property, the boost property. All stationary black hole solution is represented 
as the boost transformation of the Schwarzschild black hole. Thus various things which 
Schwarzschild black hole possess are common with stationary black holes. Using this use¬ 
ful feature we have solved the perturbation equation of the singly rotating Myers-Perry 
black hole explicitly and obtained the quasinormal mode condition, which is just an alge¬ 
braic equation, analytically. Although originally the perturbation equation is the partial 
differential equation system, the large radial gradient and boost property of the large D 
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Myers-Perry black hole reduce the equation to the analytically solvable ordinarily differ¬ 
ential equation. 

We can consider many directions of the extension of our work. In this paper we 
consider the solution with only one angular momentum in vacuum. Thus it is interesting 
to include more angular momentum, matter fields such as a gauge field or some compact 
dimensions into the effective theory of stationary black holes. These inclusions would give 
further additional degree of freedom for the deformation of the horizon. So its dynamics 
becomes much richer. Another direction of the extension is the more detail investigation 
of the effective theory of stationary black hole by searching other solution such as black 
ring solutions and going to more higher order nonlinear solutions. By constructing such 
solution we can draw the phase diagram of black holes analytically and compare with 
numerical results. Actually the higher order investigation of the effective theory for the 
static solution [9] allows very interesting phenomena such as the existence of the critical 
dimension of the non-uniform black string [36] . 
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A Singly rotating AdS Myers-Perry black hole 

In this appendix we give the metric of the AdS Myers-Perry black hole and its large D 
limit in our notation. The D = n + 3 dimensional singly rotating AdS Myers-Perry black 
hole is given by [15, 16] 13 

ds 2 = —A(r, 0) 2 dt 2 + B(r, 9) 2 dR 2 + F(r , 9) 2 {d(j) — w(r, 6)dt ) 2 

+ G(r, OfdO 2 + r 2 cos 2 9dtt 2 D _ 4 , (A.l) 


where 

2 _ A k A e (r 2 + L 2 )(r 2 + a 2 cos 2 9) 

L 2 HA^(r 2 + a 2 cos 2 9) + r 2 L 2 Ag(r 2 + a 2 ) ’ 


(A.2) 


B(r,9) 2 


r 2 R(r 2 + a 2 cos 2 9) 
n 2 R 2 


13 We change the 4 > coordinate by 


(A.3) 



from [15, 16] to satisfy the boundary condition (2.13). 
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F(r, 9) 2 = sin 2 9 


(A.4) 


r 2 + a 2 


+ 


a 2 r 2 sin 2 6 


E 2 (r 2 + a 2 cos 2 9) R 


w(r, 6) = 


ar 2 Ag 


L 2 EAAr 2 + a 2 cos 2 0) + r 2 L 2 Ag(r 2 + a 2 ) ’ 


and 


Here we defined 


and following functions 


gm ) 2 = 


r 2 + a 2 cos 2 0 

Afl 


R = 






L 2 ' 


A^ = r 2 - (r 2 + a 2 ) ^1 + j R, A# = 1 - cos 2 9. 

The horizon position, r = r+, is dehned by the positive real root of A^ by 

A k (r+) = 0. 

L is the AdS curvature and the metric is the solution of the Einstein equation 
cosmological constant 


A = - 


L 2 


At the large D limit the metric is simplified to 14 
ds 2 = 


2 dR 2 L 2 (l + a 2 cos 2 9) 


n 2 (l + a 2 )(l + L 2 )R(R — 1) 

- 1i - C ° sh2 /Wj d ? + (i + dtf 


R 


R 


2 sinhcr(0) cosh cr(9) >7 1 + a 2 cos 2 9 


R 


-dt d(j) + 


d9 z + cos 2 9dQ 2 D _ A 


As 


where we introduced 


D (l + a 2 )(l + L 2 ) 6 2 (a . 

R = ----R, cosh z a{9) = 


L 2 


(1 + a 2 )Ag 
(1 + a 2 cos 2 9)E 


dt and d<f> are defined by 


jX2 (1+^ 2 )A,9 2 j!2 1 + « 2 . 2 Q ]x2 

dt = -—- dt , dq> = ——— sm 9dcp . 


In the following we set ro = 1. 


(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 
with a 

(A.10) 

(A.ll) 

(A.12) 

(A.13) 
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At the large D limit the horizon is located on 


R 


r=r_|_ 


= 1 + 0(l/n). 


This solution can be reproduced by the leading order solution in section 2 by 

2 (1 + T 2 )A0 


Vo {oy = 


L 2 l 


lZ}±(z) = cos 9, 


and 


dz /1 + a 2 cos 2 6 2 

Te = \ - ~Ka -* R = cosh p ' 


(A. 14) 


(A.15) 


(A.16) 


The surface gravity and horizon angular velocity at the leading order of large D limit are 

(A.17) 


„ 11 + L 2 „ o 1 + L 2 

K = - -, ilH = 


2 L 2 ’ 1 + a 2 L 2 ' 

The metric (A. 11) becomes the large D metric of the Schwarzschild black hole [3] 

ds 2 = —dp' 2 — tanh 2 p dt 2 + dz 2 + sin 2 ddcf) 2 + cos 2 9dQ 2 D _ 4 , (A. 18) 

when a = 0 and l/L = 0. We changed the coordinate by R = cosh 2 p. 

Structure in 1/D correction Let us consider the 1/D correction in gu of the AdS 
Myers-Perry black hole at a = 0. Using 


^ = 1 + logR +0 (l/^ 2 ) ( A1Q ) 

n 

we find that gtt becomes 

9a = ~A(r, 6) 2 

= i+ iH + ^ +o(i/ " 2) ' (A - 20) 

At large R gtt has a linear term in p since the relation between p and R is 

R = y ^ L 2 cosh2 p- ( A - 21 ) 

Thus, in general, the higher order corrections in A(p, 6) cannot be assumed to be damping 
exponentially at the overlap region n>/)> 1 in the presence of the cosmological constant 
as mentioned in section 2. 
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A.l Ellipsoidal embedding 


The AdS-Myers-Perry black hole is described as the solution of the effective equation 
(2.39) for the ellipsoidal embedding in the AdS background. The ellipsoidal coordinate in 
AdS is given by 15 

-1(Jl i _2 _„„2 , 


d , 2 = _M^+ r V + 


L * 2 E 


(r 2 + a 2 )(L 2 + r 2 ) 

„2 , „ 2 N 2 , 


r 2 + a 2 cos 2 0 9 (r 2 + a 2 ) sin 2 0 

H- 7 -d6P + --^-( 


Afl 


The ellipsoidal embedding is defined by 


r = r( 0 ) = 


n H (e) 

cos 8 


b 2 + r 2 cos 2 9d^l 2 D _ 4 . (A. 22 ) 


(A.23) 


and the leading order metric of the AdS Myers-Perry black hole in the 1/D expansion 
is reproduced by r( 0 ) = 1 when we fix the horizon position by r + = 1 with the surface 
gravity and horizon angular velocity given in eq. (A. 17). We consider the perturbation 
around this AdS Myers-Perry black hole embedding. As done in section 2 the obtained 
perturbed solution would give the bumpy black hole solution in AdS in the perturbative 
manner. At first we fix the surface gravity to 

1 + L 2 


K = 


2 L 2 


(A.24) 


by using the time coordinate normalization in near zone. The perturbed embedding is 


77 h ( 0) = cos e\l + el1(9) + O (e 2 ) 


where 


e = Q h — 


1 + L 2 


1 + a 2 L 2 ' 

Then, perturbing eq. (2.39) with respect to e, we find the perturbative solution 


77(0) = 


a(l + a 2 ) 2 sin 2 0 


+ A 


(sin0) 5 MPAj~ 5Mp/2 


where 


£(1 — a 2 )(l + a 2 cos 2 0 ) 1 + a 2 cos 2 0 

(1 + a 2 )(r 2 - 1 ) 


(Imp = 


L 2 E 


This perturbation solution becomes regular at 0 = 0 when a = a c given by 

2 _ {/Ik — 1)L 2 + 1 

a ° ~ L 2 + 2k - 1 


(A.25) 


(A.26) 


(A.27) 


(A.28) 


(A.29) 


15 The coordinate transformation from the usual AdS coordinate to the ellipsoidal coordinate is [15] 

2 • 2 /) i-1 —1/2. 2\ • 2 /o 2 2 n 2 2 n 

r sm t) —> (r +a ) sm 0, r cos 6 1 —>■ r cos 6. 
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where A; is a positive integer. At L —>• oo eq. (A.29) reproduces the threshold angular 
momentum (2.56). Hence we expect that eq. (A.29) gives the threshold angular momen¬ 
tum of the AdS Myers-Perry black hole and branching angular momentum for the bumpy 
black hole in AdS. It is interesting to check this result by solving the perturbation equa¬ 
tion of the AdS Myers-Perr black hole and obtaining the QNM frequency as done for the 
Myers-Perry black hole in section 4. 

B Next-to-leading order analysis 

Here we briefly give results for next-to-leading order solutions. In particular we will see 
how the constancy condition of the surface gravity and horizon angular velocity, eq. (2.30), 
can be obtained. 

Our D = n + 3 dimensional metric ansatz is 
N 2 (n 7 ) 

ds 2 = --1— '-dp 2 — A(p, z) 2 dt 2 + F(p, z) 2 (d(j) — W(p, z)dt) 2 

n z 

+ G(p, z) 2 dz 2 + H(p, z) 2 qABdx A dx B . (B.l) 


We solve the Einstein equations for metric functions and the extrinsic curvature on p = 
constant surface defined by 
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K a h = , ,g ac d p g cb . 

b 2 N(p, z ) y py 

(B.2) 

The Einstein equations are decomposed on p = constant surface as 



-R + K 2 - K\K b a - {D l) ^ 2) = 0, 

(B.3) 


V a K a b - V b K = 0, 

(B.4) 


%K\ = -KK% + R\ + - ^V a V 6 lV. 

(B.5) 

The leading order solutions which are regular on the horizon are 



2 _ Aq{z) 2 Fq(z) 2 tanh 2 p 

F 0 (z ) 2 - A 0 (z) 2 Ct<f,(z) 2 tanh 2 p ’ 

(B.6) 


F(p, z) 2 = Fq(z) 2 - A 0 (z)C t<j> (z) 2 tanh 2 p, 

(B.7) 


w( , _ F 0 (z) 2 ( 1 - C tt (z)) + A 0 (z) 2 C tfl) (z) 2 Ctt(z) tanh 2 p 
K t <t,(z)Fo(z) 2 - A 0 (z) 2 C t( p{z) 3 tanh 2 p 

(B.8) 

and 

. 2ro(z) 2 (IZv(z) 1 \ , . 

0(p.^ = l „ ( KhW i2 )log(coshp), 

(B.9) 










To satisfy the boundary conditions at p 3> 1, C tt and C t $ are found to be 


Ctt(z) = 


1 


Ctt(z) = 


F 0 (z)y/1 - A 0 (zY' 


(B.ll) 


A 0 (z) 2 v A 0 (z) 

Then we can see that the leading order solutions satisfy the boundary condition at p 1 


as 


A{p,z) = \ + 0{e-»), W(p, z) = 0(e~ p ). 


(B.12) 


Furthermore the functions in this leading order solution should satisfy one equation coming 
from the vector constraint (B.4) given by 


d i 

-r lo g 

az 


Ao(z)C tt (z)\ d 

ra j- (1 - 0a( ,) )1 og 


Cu(z) 

C t4> {z) 


= 0. 


(B.13) 


This equation is rewritten by using eq. (B.ll) as 

log (A 0 (z)F 0 (z)y/l - A 0 (z) 2 ) = 0. (B.14) 


d . ... 1 — Aq(z) 2 d 

-log(4,r„M)+ Mz)l - 


We consider the 1 /D corrections to the leading order solutions. Especially, by solving eq. 
(B.5), we found that the 1/D correction to K z z , say B) K z z , has the following form 


(1) /\" z = 


/~<W 

zz 

sinh p 


+ 5K z z (p,z ), 


(B.15) 


where C^ z is an integration function with respect to p-integration. 5K z z (p, z) is a solution 
coming from the integration of the source term at next-to-leading order. cQ is used to 
eliminate the divergence of 0(p~ l ) in 5K Z z (p, z) at the horizon p = 0. Even though the 
0(p~ 1 ) divergence is eliminated by the integration function, ^K z z cannot satisfy the 
regularity condition at the horizon. Actually 5K Z z (p, z) has the following behavior after 
eliminating 0(p~ 1 ) divergence at the horizon 


5K z z (p,z) = - 


ro(z) logp 


(l-A 0 (z) 2 )A 0 (z) 2 F 0 (z) 2 p 


1 2 


F 0 (z)A' 0 (z) + (1 - A 2 )A 0 (z)F^z) + 0(1). 


(B.16) 


To satisfy the regularity condition on the horizon we should eliminate also this 0(p 1 log p) 
divergence. Then Aq(z) and Fq(z) should satisfy additional condition 


F 0 (z)A' 0 (z) + (1 - A 2 )A 0 (z)F'(z) = 0. 
This condition can be solved by 


Fo(z) = 


Vi-A 0 (z ) 2 

^hA 0 (z) 


(B.17) 


(B.18) 
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where we introduce an integration constant fin- Using this solution (B.18), eq. (B.14) is 
reduced to the equation only for Aq(z), and it can be solved by 

Aq{z) = 2kr 0 (z), (B.19) 


where k is an integration constant. 

Let us summarize above results. We obtained additional condition (B.17) from the 
boundary condition on 0(1/D) corrections. By combining this additional condition with 
the leading order vector constraint (B.14) we get the expressions of Ao(z ) and Fq(z) in 
terms of ro(z) and integration constants as eqs. (B.18) and (B.19). Using eq. (B.ll) we 
can see that eq. (2.30) is obtained by defining k = nk. 

C Spheroidal harmonics at large D 

We give some analysis for the spheroidal harmonics in the 1/D expansion. Let us consider 
the massless scalar field equation 


□4' = 0, 


(C.l) 


in the D = n + 3 dimensional flat spacetime in the ellipsoidal coordinate given by 

9 - dr 2 + (r 2 + a 2 cos 2 9)d9 2 


,9 ,9 r Z + a Z COS Z 0, 9 9 9 9 >^9 

ds 2 = -dt 2 + - lJ2 1 ,J1 1 


r z + a z 


+ (r 2 + a 2 ) sin 2 9dcj> 2 + r 2 cos 2 9dQ 2 D _ 4 . 


(C.2) 


a is an oblateness parameter. The scalar field can be decomposed as 


T = e~ iut e im4> ^(r)S(9)Y j . 


(C.3) 


Y j is the spherical harmonics on S D 4 with the angular momentum number j. Eq. (C.l) 
becomes equations for il’(r) and S(9) as 


r n—1 


d r 


„n —1 


d r + 


((r 2 + a 2 )u + am) 2 a 2 (j(j + n — 2)) 
r 2 + a 2 r 2 


-A 


i/(r) = 0, (C.4) 


and 


Q sin 6* cos n 1 0d e 

9 cos n_1 9 


sm 


(au sin 2 9 + m) 
sin 2 9 

_ j(j + n - 2) 
cos 2 9 


+ A 


5(0) = 0. (C.5) 


A is the separation constant. The spheroidal harmonics, 5(0), is defined by the solution 
of eq. (C.5). While there are numerical analysis [37] to find solution and eigenvalue of eq. 
(C.5), but it is hard to obtain analytic solutions of eq. (C.5). Here we apply the large D 
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expansion to solve eq. (C.5). At first we assume that the parameters in eq. (C.5) have 
following orders at large D 

uj = 0(1), a = 0(1), m = 0(1), j = 0(1), A = 0(n). (C. 6 ) 

The reason for A = 0(n) can be seen by considering a = 0 limit. At a = 0 eq. (C.5) is 
reduced to the equation for the spherical harmonics. The equation at a = 0 can be solved 
by S^ a=0 \d) 

S( a = °)( 0 ) = cos - 7 9 sinl m l 9 2 F 1 (—ks, ks + \m\ + j + n/ 2 , \m\ + 1 ; sin 2 9), (C.7) 

with the separation constant 

A( a= °) =l[i + n). (C. 8 ) 

ks is given by 



Then, at a = 0, the separation constant is 0(n). The spheroidal harmonics becomes the 
spherical harmonics at a = 0. Thus it is natural to assume A = 0(n) also at a / 0. 
The important observation is that the spherical harmonics (C.7) has the following large 
D behaviors 

S i l a=0) (9) = sm e ~ j 9cos j 9 + 0(l/n). (C.10) 

Usually the spherical harmonics S^ a ~°^ (9) is understood as the function which has i nodes. 
This t is the overtone number along 9 direction. At large D limit the situation is changed, 
and any nodes would disappear according to the large D limit solution (C.10). Instead 
the quantum number £ can be read from the behavior only around 9 = 0 as the degree of 
vanishing as 

S ( l a=0) (9) = 9 e (l + 0{9,n~ 1 )) . (C.ll) 

Although it is not still unclear why we can obtain the quantum number i only by the 
behavior around 0 = 0, not by the behavior both at 9 = 0 and 9 = 7 t/2 , such behavior is 
useful for the quasinormal mode analysis of the Myers-Perry black hole. 

Next we consider the contributions by a 7 ^ 0 to the spheroidal harmonics. Let us 
assume that the spheroidal harmonics S(9) is modified by the contributions from a and u 
as 

S(9) = Sj a=0) (9) (1 + SS(9)). (C. 12 ) 

Substituting this into eq. (C.5) and linearizing the equation with respect to 6S(9), we can 
obtain the equation for 5S(9). The obtained equation, however, does not have a and u 
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dependences in its equation at the leading order of the large D limit. This is because the 
contributions of a and oj are not leading order effects in eq. (C.5) in the 1/D expansion. 
Thus 5S(9) can be set to zero also for the spheroidal harmonics S(9) at the leading order 
in the 1/D expansion. Hence the spheroidal harmonics has the following large D form 

S(9) = sf = 0 \e) + O{l/n) 

= sii/ - - 7 9cos^9 + 0(l/n), (C.13) 

and the separation constant also becomes 

K = £n(l + 0{n~ 1 )) , (C.14) 

where £ is non-negative integer. The quantum number along 9 direction, £, can be obtained 
again only by the behavior around 9 = 0. If one considers the 1/D corrections to the 
spheroidal harmonics we can obtain the 1/D correction to the separation constant. In the 
analysis of this paper we consider only the leading order results of the quasinormal modes, 
so we do not pursue this analysis in more detail here. 


D Detail analysis on the vector type perturbation 


In this appendix we show the detail analysis of the vector type perturbation of the Myers- 
Perry- black hole. The vector type perturbation on the Myers-Perry black hole is decom¬ 
posed to 


h { y ) dx fi dx u = e ~ iuit e im * 


f(V)Sch Y <y)j e (i) e (A) + h { Py { ab e (A) e (s) 


(D.l) 


where the vielbein and the vector harmonics were defined in section 3. The decoupling 
perturbation variables are defined by 


y 1 + ?"o0 


(D.2) 


F (V) _ 
r 2 ~ 


cosh p 


\jl + r 0 (z) 2 sinh 2 p 


/ 2 (1 ' ) -sinhm 


1 - M z ) 2 AV) 
A > 


1 + ro(z) 2 sinh' p 


and 


(D.3) 


F (V) = f (v) p (V) = av) p (v) = H (v) 

1 1 ^ 3 J 3 5 L L 


(D.4) 


The decoupling perturbation variables are expanded by 


jf> = v 47 = v 

J- / J- jh / 


(k) p{V) 


(V) 


(k) F {v) 


k> 0 


fc >0 




(D.5) 
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Changing coordinate by R = cosh 2 p, the perturbation equation at each order becomes 


JLi 

OR 

A, 

dR 


4 

4 

d 


9R R(R %R 


(k) p( v ) 

4R(R — 1) 
0 ) F (V) 
4R(R - 1) 
9 W F (V) 

JL(k) F (y) 

SR 3 


(k) MV) 

°0 > 

(D.6) 

(fc) MV) 

°i i 

(D.7) 

(k) MV) 

°2 i 

(D.8) 

Ws™. 

(D.9) 


The boundary condition is the same with one of the scalar type perturbation, the in¬ 
going boundary condition on the horizon (4.4) and decoupling mode condition (4.5) at 
asymptotic infinity of near zone. 


Leading order (k = 0) The regular solution at the leading order is 


(0) F v = (0) F v = 


W< 


(o) 


(°) fY = (0) fY = o. 


(D.10) 


V /R ( R - !)’ 

The leading order solution has an undetermined function o) by the boundary condition. 
Next-to leading order {k = 1) The regular solution at the next to leading order is 


(i ) F v = W W 

° v/ R ( R - 1) 




w. 


( 1 ) 


v/R(R-l) 


+ upr, 


and 


( 0 ) F V = ( 1 ) pv, ( 0 ) pV = ( 1 ) pv, 


(D.ll) 


(D.12) 


where terms with the hat come from the integrations of source terms. The next to leading 
order solution has an undetermined function by the boundary condition again. At 
this order we obtain the non-trivial condition for W(q) as 


cos 9 sin 2 6( 1 + a 2 cos 2 d)W^(0) 

+ 


s 2 < 

—j( 1 + a 2 cos 2 9) + i cos 2 0(a 4 cos 2 9(oj — i) — a 3 m cos 2 9 
+ a 2 (cos 2 9(oj + i) + oj — 2i) + oj — am) Wq(6) = 0. 

This condition can be solved explicitly as 


7 a ■ (0) 5 V a 
w ia\ _ cos^sm 6 

W(q )(9) — 2 - 2~n~, 

w 1 + a z cos z (J 


(D.13) 


(D.14) 


where 


= “(1 + a 2 )(j — 1) — wm + *(1 + a 2 )oj. 


(D.15) 
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The harmonics condition on Wq{9) is eq. (4.18) as 


(D.16) 


and l is parameterized by the non-negative integer ky by 

t = j + \m\ + 2ky. (D.17) 


This harmonics condition gives the leading order QNM frequency for the vector type 
perturbation by 


^LO = 


am 

1 + a 2 



1-3 \ 

1 + a 2 J ' 


(D.18) 


Next-to-next-to leading order At this order we can obtain the solution and unde¬ 
termined function in the same manner. Furthermore we find the non-trivial condition on 
W(u(0). Although we omit the explicit form of it since it is not so illuminative, the solu¬ 
tion of the non-trivial condition gives the following behavior around 6 = 0 on the horizon 
R = 1 

% = (! + 0(6, 1/n 2 )) (D.19) 

where 

« = <0) « + (D.20) 

n 

8^ is given in eq. (D.15). is 

(1) C = cl 2 (—a 2 (j - l) 2 + 2 ia(j - l)m + j 2 - 2 j - m 2 + 1) 

— a(2a(j — 1) + im log (1 + a 2 )) + u(i(a 2 — 1) log (1 + a 2 ) 

+ 2 ia 2 (a 2 (j — 1) + iam + j — 1)) + (1 + a 2 ) 2 u 2 . (D.21) 


Then the harmonics condition gives 




This condition is solved by 


n 


v 


UJ v = UJ v o + ^o +0{l/n ^ 


where 


,v 


2 am 


W NLO - 


(t - 1 + a 2 (j - 1) + log (1 + a 2 )) 


_ i Uj _ i)2 + (! + q2 )(^ “ j)( £ + 3~ 2 ) - 2a 2 m 2 


+ j - 1 - 


(1 + a 2 ) 2 
1-3 , 2 (£-j) 


+ 


1 + a 2 (1 + a 2 ) 2 


(D.22) 


(D.23) 


log(l + o 2 ) . (D.24) 
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E Explicit forms of q(0) and di(9) 


The scalar type perturbation can satisfy the boundary condition on the horizon at the 
next to next to the leading order if the integration functions are solutions of following 
equations: 


and 


Ci(0)cfo(0) + czWcfay) + C 3 (e)c w (e) 

+ c/±(6)D" 0 ^(6) + c§(9)D'^(9) + cq(9)D^(9) = 0, 


di(e)D^(6) + d 2 (9)D'[ 0) (9) + d 3 (9)D{ 0) (9) 

+ d A (e)D {0) (6) + d b {e)c' m {9) + d 6 (0)c ( u(0) = 0. 


Each coefficients are following: 


d(e) = 


cos 2 9 sin 3 9 (l + a 2 cos 9 2 ) 

2(1 + a 2 ) 5/2 


(E.l) 


(E.2) 


(E.3) 


, cos 9 sin 2 9 r 9 A „ , 9 , . x 

c 2 (9) = - rj- 2 a cos 9 [a (j — iu) + iam — iui + 2) 

2(1 + a 2 ) 5/2 L V ’ 

+ cos 2 9 (a 2 (4j — 2 iuj — 5) + 2 iam — 2 ioo — 2) + 2j + 1 


(E.4) 


C3(6>) = 


sin 


’ 2 cos 6 9 ( a 2 (j — *w) + iam — + 2)' 


2(1 + a 2 ) 5/2 

+ cos 2 0 (a 2 (3j 2 + j( — 4 — 2iuj) — l) + 2 iajm + j( — 4 — 2ioj) + 4) 

+ cos 4 0^a 4 (3j 2 + j(— 6 — 4iuj) — cu 2 + + 4) + 2a 3 m(2ij + a; — 3i) 

— a 2 (m 2 + 2 (2 i(j — 2)u + w 2 — l)) + 2am(uj — i) — w(cj — 2i)^ 


+ .r + 2.7 - 3 


(E.5) 


c 4 (0) = - 


sm 


cos 2 9 (a 2 (—5j + 3zcu + 2) — 3zam + j + 3zw + 5) 


(1 + a 2 ) 4 cos 9 

+ cos 4 0 (a 4 (—2j + 3iu + 3) — 3 ia 3 m + a 2 (j + 2 iw + 1) — 
— za 2 (l + a 2 ) cj cos 6 0 — 3 j — 5 


(E.6) 
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c 5 (0) = 


sin 9 


cos 2 9 (a 2 (lOj 2 + 2j(8 — 3zo>) — 7iu — 3) 


(1 + a 2 ) 4 cos 2 9 (1 + a 2 cos 2 9) L 
+ za( 6 j + 7 )m — 2 j 2 + j(—19 — 6 za;) — 7iu — 9^ 

+ a 4 (a 2 + l) a; cos 10 9 (a 2 ( 2 zj + 2 a; — i) — 2am + 2uj + 3z) 

+ cos 4 9 (a 4 (l2j 2 + j {—6 — 16 iu) — 3a ; 2 — ioj — 9) + a 3 m(16zj + 6 a; + z) 

— a 2 ( 6 , 7 2 + 4j(7 + 3za>) + 3m 2 + 6 a > 2 — 9 iui + 2) 

+ am(—2ij + 6 a; — 9z) + j(5 + 4iu) — 3 a ; 2 + lOza; + 4^ 

+ a 2 cos 8 9 (a 6 (j 2 + j{—3 — 4 iuj) — 3a ; 2 + 4za; + 2) + 2a 5 m(2ij + 3a; — 2i) 

— a 4 (2j 2 + j{—2 — 4za>) + 3 m 2 + 2a>(a; + z)) + a 3 m(—2zj + 2a; + z) 

+ a 2 ( 8 zja; + j + 5a ; 2 — 5za; — 2) + am(—4a; + z) + a;(4a; + z)^ 

+ cos 6 9 (a 6 ( 6 j 2 + j {—12 — 14za;) — 6 a ; 2 + lOza; + 5) + 2a 5 m(7ij + 6 a; — 5z) 

— a 4 ( 6 , 7 2 + j(7 + 4za>) + 6 m 2 + 10a ; 2 — 15za; — 15) + 2a 3 m(—2ij + 5a; — 4z) 
+ a 2 (2j(3 + 5za>) — 2a ; 2 + 3za> + 6 ) + am(—2a> + i) + 2a;(a; — z)^ 


+ 3 j 2 + 13 j + 5 


di{6) = - 


cos 9 sin 3 9 


d 2 {9) = - 


(1 + a 2 ) 2 (1 + a 2 cos 2 9) 
sin 2 9 


(E.7) 


(E. 8 ) 


(1 + a 2 ) 2 (1 + a 2 cos 2 9) 2 - 


3.7 + 1 


+ cos 2 9 (a 2 ( 6.7 — 3za; — 4) + 3z am — 3iuo — l) 
+ a 2 cos 4 9 (3a 2 (j — za; — 1) + 3zam — 3iu + l) 


(E.9) 


sin 9 


cos 2 #^a 2 (l2j 2 — 6iju + j + ico + 3) 


dz{6) =--,, 

(1 + a 2 ) 2 cos# (1 + a 2 cos 2 #) 3 

+ za(6j — l)m — 6ijui — 6 j + za; — 1^ + a 4 cos 8 9 ^a 4 (3 j 2 — 6ijui — 4 j — 3a; 2 + 5za; + 2) 

+ a 3 m(6ij + 6a; — 5z) + a 2 (—6zja; + 4 j — 3m 2 — 6a; 2 + 2za> — 2) + 3am{2u + z) 

— 3a;(a; + z)^ + a 2 cos 6 #^a 4 (l2j 2 — 18zja; — 17j — 6a; 2 + 15za; + 7) + 3a 3 m{Qij + 4a; — 5z) 

+ a 2 (—18zja; + 2 j — 6m 2 — 12a; 2 + 13za; + 5) + 2am(6ui + z) — 6a> 2 — 2za> + 8^ 

+ cos 4 #^a 4 (l8j 2 + j {—17 — 18za;) — 3a; 2 + llza; — l) + a 3 m(18zj + 6a; — llz) 

— a 2 (2j(4 + 9za;) + 3m 2 + 6a; 2 — 12za; + 10) + am(6a; — z) + (—3a; + z)a>^ 


+ 3 j 2 + 5.7 + 1 


(E.10) 


dm = 


VT+ a 2 cos # 2 sin 2 9 
(1 + a 2 ) (1 + a 2 cos 2 9) ’ 


(E.ll) 
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dm = 


cos 6 sin 9 (cos 2 9 (a 2 (j — iu — 2) + iam — iu) + j — l) 


\/l + a 2 (1 + a 2 cos 2 0) 

The coefficients cg(0) and d/#) have very lengthy forms. They are written as 

k =7 


and 


where 


c&{9) = 


1 


( a 2 + l) 4 cos 3 9 (1 + a 2 cos 2 9) z f/ 


E 


2/^6 




Cfl COS 


k=0 


k=6 


dm = 


(1 + a 2 ) 2 cos 2 9 (1 + a 2 cos 2 #) 4 /" 


d\ cos 2k 9, 


k =0 


4 U) = i (J 2 + 8 j + 12) , 

Cg 1 * = j (5j 2 — 3 iju + 23 j — 1C )iu + 16) 

+ ia(3j + 10)m — j 2 — 3 ij(u — 5 i) — 10 iu — 28^ , 


(E.12) 

(E.13) 

(E.14) 

(E.15) 

(E.16) 


4 2) = a 4 (lOj 3 + j 2 (18 - 13iu) - j (3c/ + 19*a; + 8) - 2 (a; 2 - 2 iu + 4)) 

+ a 3 m (13 ij 2 + 6jui + 19ij + 4a; — 4 i) — a 2 ^5j 3 + 5/(9 + 2iu) 

+ j (3m 2 + 6c/ + 44) + 2 (m 2 + 2c/ — 2 iu + 4)^ + am (—2 ij 2 + j(6co — 17 i) + 4 u) 

+ /(7 + 3 iu) + j (-3c j 2 + 19iu + 20) - 2c/, (E.17) 


4 3) = a 6 (10/ +/(—4 - 22ia>) - j (lie/ + 14) + i (u 3 + 8a; + 8z)) 

+ a 5 m (22 ij 2 + 22 ju — i (3c/ + 8)) + a 4 ^—10/ +/(—42 — 9iu) 

+ j (—11m 2 + (—19a; + 39*)a>) + % (3m 2 + 4) u + 3iu 3 + 4a; 2 + 8^ 

— ia 3 m (8/ + 2/19 + 9 iu) + m 2 + 6a; 2 — 8iu — 4) + a 2 ^/(22 + 13iu) 

+ j (m 2 — 5a; 2 + 30ia; + 34) + m 2 (4 + 37a;) + 3iu 3 + 8a; 2 — 4za; + 16^) 

— am(/ 4a; — 6 i) + (8 + 3iu)u) + 3ju 2 — 9iju — 4 j + iu 3 + 4a; 2 , (E.18) 


4 4) = a 8 (5/ + /(—10 — 18ia;) — 3ju(5u — 6i) + 3iu 3 + 7a; 2 + 8) 

+ a 7 m (18*/ + 6j(5a; — 3 i) + (—14 — 9iu)u) + a 6 ^—10/ + /(—6 + 47a;) 

+ j (—15m 2 — 19c/ + 28ia; + 18) + m 2 (7 + 97a;) + 8ic/ + 15a; 2 — 10ia; + 32^ 

— ia 5 m (12/ + 4/5 + 4iw) + 3m 2 + 16c/ - 21ia- - 12) + a 4 (/(24 + 227a;) 

+ j (3m 2 + 7c / — 107a; + 8) + m 2 (6 + 8kS) + 6iu 3 + 7 c/ — 12zca + 16^ 

— a 3 m(14j(a; — i) + (4 + 5iu)u) + a 2 (j (lla; 2 — 20ia> — 6) + m 2 (—1 — iu) 

— 3c/ — 2iu — 8^ + am(3 + 2iu)u + (—2 — iu)u 2 , (E.19) 
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Cg 5) = a 2 |^a 8 (j 3 + j 2 (— 3 — 7ia) + j (—9a 2 + 10ia + 2) + a (3ia 2 + 6a — 4i)) 

+ a 7 m (7 ij 2 + 2j(9a - 5 i) - i(-3u + 2 i) 2 ) + a 6 (-5 j 3 + j 2 { 9 + Ilia;) 

— 3j (3m 2 + a(a + 4z)) + m 2 (6 + 9ia) + 6ia 3 + a; 2 — 4ia — 8^ 

— ia 5 m (8 j 2 — 6 j + 3 m 2 + 12a; 2 — ia — 4) + a 4 ^2j 2 (5 + 9ia) 

+ j (3m 2 + 21a; 2 - 32ia - 6) + 2i (3m 2 - l) a - 17a; 2 - 24^ 

+ a 3 m (—2j(9a — 5i) + 3ia 2 + 14a; — 4i) + a 2 ^5j(3a — 2i)a + m 2 (—2 — 3ia) 

— 6ia 3 — 13a; 2 — 16^ + 3am(l + 2ia)a + a (—3ia 2 — a; + 2i) j, (E.20) 


Cg 6 ^ = a 4 ^a 8 a (— ij 2 — 2ja + ij + iu 2 + a;) + a 7 m ( ij 2 + j(4a — i) + (—2 — 3ia)a) 

+ a 6 (— j 3 + j 2 (3 + 6ia) — j (2m 2 — 5a; 2 + 12ia + 2) + m 2 (l + 3ia) + (—7a; + 6i)a) 

— ia 5 m (2 j 2 + j (—5 — 6ia) + m 2 + 9ia + 4) + a 4 (j' 2 (l + 7ia) + j (m 2 + a;(16a; — 9*)) 

— 2m 2 + a; (—6ia 2 — 13a; + 4i)^ + a 3 m (j(—10a + 2i) + 9ia 2 + 8a — 4i) 

+ a 2 (j (9a 2 + 4ia + 2) + m 2 (—1 — 3ia) + a (—8ia 2 — a + 2i)) 

+ 3iama(2a + i) + a (—3ia 2 + 4a + 4i) J, (E.21) 


4 7) = a 6 (a 2 + l) a^a 4 (ij 2 + j(2a — i) + (—1 — ia)a) + a 3 m(—2j + 2ia + 1) 

+ a 2 (j(2a + 3i) — i (m 2 + 2a 2 + 2ia + 2)) + am(—3 + 2ia) 

-ia 2 + 3a + 2i), (E.22) 


4 0) = ~j(3 + 2) 2 , 


(E.23) 


44 = j (a 2 (—6 j 2 + j(—11 + 3ia) + 2ia — 8 ) — ia(3j + 2)m 


(E.24) 


df ] = a 4 (—15j 3 +j 2 (- 3 + 15ia) + j (3a 2 - 5ia + 2) - 2 (a 2 - 2ia + 4)) 

+ a 3 m (—15ij 2 — 6ja + 5 ij + 4a — 4i) + a 2 (j 2 (17 + 15ia) + j (3m 2 + 6a 2 — lOia + 26) 

— 2 (m 2 + 2a 2 — 2ia + 4)^ + am(4a + j(— 6a + 5i)) + j (3a 2 — 5ia — 4) — 2a 2 , (E.25) 


44 


= a 6 (—20 j 3 + 6j 2 (3 + 5 ico) + 2 j (6a; 2 — 15 ioj + 5) — za; (a; 2 — lOza; — 4)) 

+ a 5 m (—30zj 2 — 6j(4a; — 5i) + 3za> 2 + 20a; — 4z) + a 4 ^6j 2 (3 + 5za>) 

+ (3m 2 + 6a; 2 — 10za> + 6) + m 2 (—10 — 3iuj) — 3za> 3 — 20a; 2 — 8za; + 16^ 

+ a 3 m (j(—24a; + 10?) + z (m 2 + 6a; 2 — 20za> + 12)) + a 2 ^2j (6a; 2 — 5za; — 9) 

— 3z (m 2 + 4) u — 3za; 3 — 10a; 2 + 16^ + 3zama; 2 — ico 3 , (E.26) 


df ] = -a 2 j^a 6 (l5j 3 + j 2 (-22 - 30ia;) - 2 j (9a ; 2 - 20za; + l) + 3ia ; 3 + 16a; 2 - 4ia; + 8 ) 

+ a 5 m (30zj 2 + 4j(9a; — lOz) — 9za ; 2 — 32a; + 4z) + a 4 (j' 2 (—2 — 30za;) 

— 2 j (9m 2 + 18a; 2 — 20za; + 2) + m 2 (16 + 9za;) + 9za ? 3 + 30a> 2 + 8 za; + 16^ 

— a 3 m (—36ja; + 3 im 2 + 18za; 2 + 28a; + 12z) + a 2 (j (26 — 18a; 2 ) + m 2 (—2 + 9za>) 

+ 9za ; 3 + 12a ; 2 + 4za; + 16^ + am (—9za ; 2 + 4a; + 8 z) + 3za ; 3 — 2a ; 2 — 8iui + 8 ^ , (E.27) 


(I 4 = a 4 ^a 6 (— 6 j 3 + 3j 2 (3 + 5zaz) + 2 j ( 6 a > 2 — lOza; — l) + u (—3za ; 2 — 10a; + 4z)) 

+ a 5 m (—15 ij 2 — 4j(6a; — 5z) + 9za ; 2 + 20a; — 4z) + a 4 ^j 2 (—7 + 15za;) 

+ 2 j ( 6 m 2 + a;(12a; — 5z)) + m 2 (—10 — 9za>) — 9iu 3 — 16a; 2 + 8 ^ 

+ a 3 m (—2j(12a; + 5z) + 3im 2 + 18za> 2 + 12a; + 4z) + a 2 ^2j ( 6 a > 2 + 5za; — 7) 

+ m 2 (4 — 9za>) — 9za ; 3 — 2a > 2 + 4za; + 16^ + iam (9a ; 2 + 8iu — 8 ) 

— 3za ; 3 + 4a ; 2 + 8iuo + 8 j , (E.28) 


d^ = —a 6 ^a 6 (j 3 + j 2 (—1 — 3zaz) — 3ja;(a; — z) + (2 + iuj)uj 2 ) 

+ a 5 m (3zj 2 + j(6a; — 3z) + (—4 — 3za;)a;) + a 4 ^j 2 (3 — 3za>) 

— j (3m 2 + 6a> 2 + 2za; + 2) + m 2 (2 + 3za>) + a; (3za; 2 + 2a> + 4z)^ 

— ia 3 m {j (—5 + 6za>) + m 2 + 6a> 2 + 4) + a 2 (j (—3a; 2 — 5za> + 2) 

+ z (m 2 (3a> + 2z) + a; (3a; 2 + 2za> + 4))^ + am(4 — 3za>)a; + za; 2 (a; + 2z) j. (E.29) 
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F Explicit form of 5% 

The QNM condition for the scalar type perturbation of the Myers-Perry black hole is 


st=e-j. 


This condition can be rewritten as 




where 


G(lo) = (l + a 2 ) 3 a 4 (j — iu) 2 + 2a 3 m(uj + ij ) 

— a 2 (j(—2£ + 2iuo + 4) + 2i£uj + m 2 + 2co 2 — Aiu> — 4) 
+ 2 am(i7 + u — 2i) + (7 — iu — 2) 2 


and 


Each coefficient is 


-F(w) = f 5 u 5 + / 4 w 4 + / 3 w 3 + f 2 u 2 + fiu + / 0 . 

/ 5 = i(l + a 2 ) 6 , 

/ 4 = — (l + a 2 ) 1 (a 2 (5j — 2) + 5mm + 57 — 6) , 


(F.l) 


(F.2) 


(F.3) 

(F.4) 

(F.5) 

(F.6) 


/ 3 = -i (1 + a 2 ) 3 (a 4 j(10j - 7) + 4m 3 (5j - 2)m + a 2 (j(207 - 21) 

- 97 - 10m 2 + 12) + 2m(107 - ll)m + 107 2 - 237 + 12^ , (F.7) 


h = (1 + a 2 ) 2 (o 6 j 2 (10j - 9) + 3m 5 j(10j - 7)m + a 4 (3j 2 (107 - 9) 

- 6 j (47 + 5m 2 - 5) + 2 (7 + 6 m 2 - 4)) - ia 3 m(j( 57 - 607) + 277 
+ 10m 2 - 30) + a 2 (30 j(£ - l) 2 - 157 2 - 307m 2 + 347 + 30m 2 - 16) 

+ 3m (107 2 - 217 + 10) m + 107 3 - 337 2 + 327 - 8^, (F.8) 


/i = i (1 + a 2 ) (5a 8 (j - l)j 3 + 2ia 7 j 2 (10j - 9 )m + a 6 j(5j 2 (47 - 3) 

- 3 j (77 + 10m 2 - 8) + 47 + 21m 2 - 12) - 4m 5 m(-3j 2 (57 - 4) 

+ j (127 + 5m 2 - 12) - 7 - 2m 2 + 3) + a 4 (3j 2 (107 2 - 177 + 8) 

+ j (—277 2 + 7 (56 - 60m 2 ) + 51m 2 - 24) + 47 2 + 37 (9m 2 - 4) + m 2 (5m 2 - 24)) 

- 2 ia 3 m (-6 j (57 2 - 97 + 4) + 157 2 + 27 (5m 2 - 14) - 9m 2 + 12) 

+ a 2 (j (207 3 - 577 2 + 527 - 12) - 117 3 + 7 2 (32 - 30m 2 ) + 37 (19m 2 - 8) - 24m 2 ) 

+ 4m (57 3 - 157 2 + 137 - 3) m + 7 (57 3 - 217 2 + 287 - 12)) , (F.9) 
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and 


/o = - ( a 2 j + iam + 7) 2 (a 6 (j - 1 )j 2 + 3 ia 5 (j - l)jrn 
+ a 4 (3j 2 (7 - 1) + j (-47 - 3m 2 + 6) + 2 (7 + m 2 - 2)) 

- ia 3 m (j(7 - 67) + 57 + m 2 - 6) + a 2 (j (37 2 - 87 + 6) - 37 2 

- 37m 2 + 107 + 4m 2 - 8) + 3za (7 2 - 37 + 2) m + (7 - 2) 2 (7 - 1)) . (F.10) 


References 

[1] V. Asnin, D. Gorbonos, S. Hadar, B. Kol, M. Levi and U. Miyamoto, “High and Low 
Dimensions in The Black Hole Negative Mode,” Class. Quant. Grav. 24 (2007) 5527 
[arXiv:0706.1555 [hep-th]]. 

[2] R. Emparan, R. Suzuki and K. Tanabe, “The large D limit of General Relativity,” 
JHEP 1306 (2013) 009 [arXiv: 1302.6382 [hep-th]]. 

[3] R. Emparan, D. Grumiller and K. Tanabe, “Large D gravity and low D strings,” 
Phys. Rev. Lett. 110 (2013) 251102 [arXiv: 1303.1995 [hep-th]]. 

[4] R. Emparan and K. Tanabe, “Holographic superconductivity in the large D expan¬ 
sion,” JHEP 1401, 145 (2014) [arXiv:1312.1108 [hep-th]]. 

[5] R. Emparan and K. Tanabe, “Universal quasinormal modes of large D black holes,” 
Phys. Rev. D 89 (2014) 064028 [arXiv:1401.1957 [hep-th]]. 

[6] R. Emparan, R. Suzuki and K. Tanabe, “Instability of rotating black holes: large D 
analysis,” JHEP 1406 (2014) 106 [arXiv:1402.6215 [hep-th]]. 

[7] R. Emparan, R. Suzuki and K. Tanabe, “Decoupling and non-decoupling dynamics 
of large D black holes,” JHEP 1407 (2014) 113 [arXiv: 1406.1258 [hep-th]]. 

[8] R. Emparan, R. Suzuki and K. Tanabe, “Quasinormal modes of (Anti-)de Sitter black 
holes in the 1/D expansion,” JHEP 1504, 085 (2015) [arXiv: 1502.02820 [hep-th]]. 

[9] R. Emparan, T. Shiromizu, R. Suzuki, K. Tanabe and T. Tanaka, “Effective theory 
of Black Holes in the 1/D expansion,” JHEP 1506, 159 (2015) [arXiv: 1504.06489 
[hep-th]]. 

[10] S. Bhattacharyya, A. De, S. Minwalla, R. Mohan and A. Saha, “A membrane 
paradigm at large D,” arXiv: 1504.06613 [hep-th], 

[11] T. Damour, “Surface Effects in Black Hole Physics”, Proceedings of the Second Marcel 
Grossmann Meeting on General Relativity, (edited by R. Ruffni, North Holland, 1982) 
p. 587. 


47 


[12] R. H. Price and K. S. Thorne, “Membrane Viewpoint on Black Holes: Properties and 
Evolution of the Stretched Horizon,” Phys. Rev. D 33 (1986) 915. 

[13] R. Emparan, T. Harmark, V. Niarchos and N. A. Obers, “Essentials of Blackfold 
Dynamics,” JHEP 1003 (2010) 063 [arXiv:0910.1601 [hep-th]]. 

[14] R. C. Myers and M. J. Perry, “Black Holes in Higher Dimensional Space-Times,” 
Annals Phys. 172, 304 (1986). 

[15] S. W. Hawking, C. J. Hunter and M. Taylor, “Rotation and the AdS / CFT corre¬ 
spondence,” Phys. Rev. D 59, 064005 (1999) [hep-th/9811056]. 

[16] G. W. Gibbons, H. Lu, D. N. Page and C. N. Pope, “The General Kerr-de Sitter 
metrics in all dimensions,” J. Geom. Phys. 53, 49 (2005) [hep-th/0404008]. 

[17] O. J. C. Dias, J. E. Santos and B. Way, “Rings, Ripples, and Rotation: Connecting 
Black Holes to Black Rings,” JHEP 1407, 045 (2014) [arXiv: 1402.6345 [hep-th]]. 

[18] R. Emparan, P. Figueras and M. Martinez, “Bumpy black holes,” JHEP 1412, 072 
(2014) [arXiv: 1410.4764 [hep-th]]. 

[19] B. Kleihaus, J. Kunz and E. Radu, “Black rings in six dimensions,” Phys. Lett. B 
718, 1073 (2013) [arXiv: 1205.5437 [hep-th]]. 

[20] J. M. Bardeen, B. Carter and S. W. Hawking, “The Four laws of black hole mechan¬ 
ics,” Commun. Math. Phys. 31, 161 (1973). 

[21] M. M. Caldarelli, O. J. C. Dias, R. Emparan and D. Klemm, “Black Holes as Lumps 
of Fluid,” JHEP 0904, 024 (2009) [arXiv:0811.2381 [hep-th]]. 

[22] W. Hsiang, Z. Teng, W. Yu, “Examples of constant mean curvature immersions of 
the 3-sphere into euclidean 4-space,” Proc. Natl. Acad. Sci. USA, 79 (1982) 3931. 

[23] G. W. Gibbons, D. Ida and T. Shiromizu, “Uniqueness and nonuniqueness of 
static black holes in higher dimensions,” Phys. Rev. Lett. 89 (2002) 041101 [hep- 
th/0206049], 

[24] R. Emparan and R. C. Myers, “Instability of ultra-spinning black holes,” JHEP 0309, 
025 (2003) [hep-th/0308056]. 

[25] O. J. C. Dias, P. Figueras, R. Monteiro, J. E. Santos and R. Emparan, “Instability 
and new phases of higher-dimensional rotating black holes,” Phys. Rev. D 80, 111701 
(2009) [arXiv:0907.2248 [hep-th]]. 

[26] O. J. C. Dias, P. Figueras, R. Monteiro and J. E. Santos, “Ultraspinning instability 
of rotating black holes,” Phys. Rev. D 82, 104025 (2010) [arXiv: 1006.1904 [hep-th]]. 


48 


[27] H. Kodama and A. Ishibashi, “A Master equation for gravitational perturbations of 
maximally symmetric black holes in higher dimensions,” Prog. Theor. Phys. 110 , 701 
(2003) [hep-th/0305147]. 

[28] D. Gorbonos and B. Kol, “A Dialogue of multipoles: Matched asymptotic expansion 
for caged black holes,” JHEP 0406 , 053 (2004) [hep-th/0406002], 

[29] O. J. C. Dias, G. S. Hartnett and J. E. Santos, “Quasinormal modes of asymptotically 
flat rotating black holes,” Class. Quant. Grav. 31 (2014) 24, 245011 [arXiv:1402.7047 
[hep-th]]. 

[30] H. Kodama and A. Ishibashi, “Master equations for perturbations of generalized 
static black holes with charge in higher dimensions,” Prog. Theor. Phys. Ill, 29 
(2004) [hep-th/0308128]. 

[31] H. Kodama, R. A. Konoplya and A. Zhidenko, “Gravitational stability of simply 
rotating Myers-Perry black holes: Tensorial perturbations,” Phys. Rev. D 81 , 044007 
(2010) [arXiv:0904.2154 [gr-qc]]. 

[32] R. Gregory and R. Laflamme, “Black strings and p-branes are unstable,” Phys. Rev. 
Lett. 70 , 2837 (1993) [hep-th/9301052], 

[33] H. Kudoh, “Origin of black string instability,” Phys. Rev. D 73 , 104034 (2006) [hep- 
th/0602001], 

[34] G. S. Hartnett and J. E. Santos, “Non-Axisymmetric Instability of Rotating Black 
Holes in Higher Dimensions,” Phys. Rev. D 88, 041505 (2013) [arXiv: 1306.4318 [gr- 
qc]]. 

[35] M. Shibata and H. Yoshino, “Bar-mode instability of rapidly spinning black hole in 
higher dimensions: Numerical simulation in general relativity,” Phys. Rev. D 81 , 
104035 (2010) [arXiv: 1004.4970 [gr-qc]]. 

[36] R. Suzuki and K. Tanabe, “Non-uniform black strings and the critical dimension in 
the 1/D expansion,” arXiv:1506.01890 [hep-th], 

[37] E. Berti, V. Cardoso and M. Casals, “Eigenvalues and eigenfunctions of spin-weighted 
spheroidal harmonics in four and higher dimensions,” Phys. Rev. D 73, 024013 (2006) 
[Phys. Rev. D 73, 109902 (2006)] [gr-qc/0511111]. 


49 


